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Synopsis 


Insertion of a vibration isolator between the source and receiver of vibration is a very 
common method of vibration control. The function of an isolator is to reduce the 
magnitude of motion transmitted from the foundation to the receiver or to reduce the 
magnitude of force transmitted from the source to the foundation. 

Shock refers to the transient response following a sudden change in some motion 
variable as experienced by, for example, an aircraft while landing, a vehicle while 
encountering a bump on the road or the foundation of a forging machine when the 
hammer is dropped on the anvil. The simplest idealization of shock isolation system 
is a single degree-of-freedom (SDOF) model in which the masses are considered to be 
concentrated and the flexible members are considered to be massless. Then, the isolator 
is assumed to be massless with linear stiffness and linear damping. The foundation is 
taken as ideally rigid. Analysis of such linear shock isolation systems have been widely 
reported in the literature. 

A large class of isolator systems is also subjected to random excitation. For example, 
the landing gear of an aircraft (while taxiing on the runway) and a vehicle suspension 
are subjected to track induced random excitations. Isolators of ground base structures 
are subjected to earthquakes which again can be modelled as random excitations. In 
the most general case, the excitation can be a combination of transient, periodic or 
random disturbances. A brief review of previous work on shock isolators and random 



vibration isolators especially with reference to vehicle vibration control is provided in 
the first chapter of the thesis. 

Nonlinearity is ubiquitous in nature. Linearity is an approximation to reality. In 
shock and vibration systems, isolators such as air springs, elastomeric dampers and 
wire-rope isolators are inherently nonlinear. Assumption of Hooke’s law for springs 
and linear viscous damping for dampers is done just for mathematical simplicity. 
Sometimes, the amplitudes of steady-state vibration are small enough to justify the 
assumption of linearity. However, the transient displacements and velocities may often 
be sufficiently large when the nonlinearity in springs and dampers cannot be ignored. 

The primary objective of this thesis is to analyze the role of nonlinearity in both 
the stiffness and damping on the performance of shock isolation systems. 

In chapter 2, base excited nonlinear SDOF isolator model has been considered. 
A simple technique combining the straightforward perturbation method with Laplace 
transform is used for obtaining solution in closed form. This method is applicable 
for any order of nonlinearity, both in the restoring and damping forces, expressed 
in the form of polynomials. However, the results are obtained with only cubic non- 
linearity both in the restoring and damping forces. Three types of base excitations, 
namely, the rounded step, the rounded pulse and the oscillatory step are considered. 
Analytical results obtained by the present method are compared with those obtained 
by direct numerical integration. It has been shown that a critical value of the positive 
coefficient of the cubic damping term maximizes the peak values of relative and absolute 
displacements. This is true for any power-law damping force with an index greater 
than 1. Numerical results are included for a typical elastomeric damper which can 
be modelled by a combination of linear viscous damping and a dissipative quadratic 
damping with a negative coefficient. It is observed that the overall performance of 
a shock isolator improves if the nonlinear damping term is symmetric and quadratic 
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with a negative coefficient as exhibited by elastomers up to a certain value of the 
amplitude of deformation. It has been established that the nonlinearity in damping 
predominantly affects the shock isolator response, whereas the effect of nonlinearity in 
stiffness is insignificant. So, in the subsequent chapters only the damping (and not the 
stiffness) is taken to be nonlinear. 

The performance of the isolator system considered in chapter 2 is re-examined 
under force shock excitations. These results are included in chapter 3. Unlike with 
base excitation, the presence of nonlinearity in the damping force is seen to marginally 
improve the performance of a linear isolator if the system is excited by a shock force 
on the mass. 

In the presence of inherent nonlinear cubic damping in a nonlinear shock isolator 
with base excitation, four different modifications of the isolator are considered to im- 
prove its performance. These modifications are (i) an isolator with a Coulomb damper, 
(ii) a three element isolator with the nonlinear damper mounted elastically, (iii) an iso- 
lator with vibration absorber and (iv) a two stage isolator. Overall, it is seen that the 
three element and two stage isolators are preferable in the presence of nonlinear cubic 
damping. These results are presented in chapter 4. 

In chapter 5, a two degree-of-freedom (2-DOF) model for a vehicle undergoing 
heave and pitch motions is investigated for shock excitations. It is assumed that the 
vehicle is moving with a constant speed so that the rear wheels are subjected to the 
same input as the front wheels but with a time delay. Both coupled and decoupled 
symmetric systems are considered. It is shown that the decoupled symmetric system 
gives optimum isolation performance. Like in SDOF systems, nonlinear damping is 
detrimental to 2-DOF systems for shock excitation with a two element isolator. A 
three element isolator nullifies the adverse effect of the nonlinear damping in the 2- 
DOF system as well. 
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The vehicle model considered in chapter 5 is investigated again with random exci- 
tation. This study is presented in chapter 6. Here, the front and rear suspensions are 
subjected to the same filtered white noise but with some time delay governed by the 
speed of the vehicle. Shinozuka’s spectral method is used for generation of the random 
process. Monte-Carlo simulation has been employed to find the response of the vehicle. 

The major conclusions of the present thesis are listed below. 

(i) A simple analytical method has been developed using perturbation technique in 
conjunction with Laplace transform to determine the transient response of an SDOF 
system in the presence of nonlinear, dissipative shock isolators. It is observed that 
in most cases, the closed form perturbed solution matched well with the numerically 
integrated results. 

(ii) It has been established that the effect of nonlinear damping is more predominant 
than the effect of nonlinear stiffness so far as the performance of shock isolators for 
base excitation is concerned. The presence of nonlinear n th power-law damping, with 
a positive coefficient and the index n greater than unity, is harmful for base excited 
shock isolators. However, quadratic damping with a negative coefficient improves the 
performance of a shock isolator. 

(iii) With force excitation, the presence of nonlinearity in damping does not deterio- 
rate the performance of a shock isolator, rather it marginally improves the performance 
of a linear shock isolator. 

(iv) To improve the performance of a shock isolator with base excitation, in the 
presence of nonlinear cubic damping, four methods, namely, Coulomb damped isolator, 
three element isolator, vibration absorber and two stage isolator have been considered. 
Overall, a three element isolator and a two stage isolator are preferable in the presence 
of nonlinear cubic damping. 

(v) 2-DOF coupled nonlinear asymmetric systems behave like a symmetric system 



for a particular choice of parameter values. At very high speeds and for forward centre 
of mass, high damping and high stiffness at the front suspension rather than at the 
rear suspension give the optimum response for the decoupled system. 

(vi) The detrimental effect of nonlinear cubic damping in a two element isolator 
and its mitigation by using a three element isolator are exhibited by both SDOF and 
2-DOF models. 

(vii) With a random base excitation, ‘Stochastic Resonance’ type phenomenon is 
observed for both heave and pitch displacements when the damping is linear. Unlike 
for shock excitations, in this case the nonlinear cubic damping reduces heave and pitch 
displacements. However, the presence of nonlinearity in the damping force increases 
both heave and pitch accelerations. 



Contents 


Certificate i 

Acknowledgements ii 

Synopsis iv 

Table of Contents ix 

List of Figures xiii 

List of Symbols xviii 

1 INTRODUCTION 1 

1.1 Introduction 1 

1.2 Review of Previous Work '4 

1.2.1 Shock Isolators 4 

1.2.2 Random Vibration Isolators 6 

1.3 Objectives and Organization of the Present Work 11 

2 RESPONSE OF NONLINEAR DISSIPATIVE SHOCK ISOLATORS 

TO DISPLACEMENT EXCITATION 14 

2.1 Introduction . 14 

2.2 Equations of Motion 17 

2.3 Types of Input Displacements of The Base 17 


ix 



2.4 Closed Form Solutions 24 

2.4.1 Zeroth Order Solution 25 

2.4.2 First Order Solution 29 

2.4.3 Second Order Solution 31 

2.5 Results and Discussions 31 

2.5.1 Validation of Closed Form Solutions 31 

2.5.2 Effect of Nonlinear Damping 38 

3 RESPONSE OF NONLINEAR DISSIPATIVE SHOCK ISOLATORS 

TO FORCE EXCITATION 45 

3.1 Introduction 45 

3.2 Equations of Motion 46 

3.3 Results and Discussions 48 

3.3.1 Effect of Nonlinear Damping 51 

4 PERFORMANCE OF NONLINEAR ISOLATORS AND 

ABSORBERS TO SHOCK EXCITATIONS 58 

4.1 Introduction 58 

4.2 Equations of Motion 60 

4.2.1 Type 1 Isolator : Shock Isolator with a Coulomb Damper .... 60 

4.2.2 Type 2 Isolator : Three Element Shock Isolator 60 

4.2.3 Type 3 Isolator : Transient Vibration Absorber 63 

4.2.4 Type 4 Isolator : Two Stage Shock Isolator 65 

4.3 Results and Discussions 67 

4.3.1 Type 1 Isolator 70 

4.3.2 Type 2 Isolator 71 

4.3.3 Type 3 Isolator 75 

4.3.4 Type 4 Isolator 81 


x 



5 2-DOF HEAVE-PITCH MODEL OF A VEHICLE FOR SHOCK 

EXCITATION 86 

5.1 Introduction 86 

5.2 System Models 87 

5.2.1 Model I 87 

5.2.2 Model II 91 

5.3 Theoretical Analysis for Special Cases 93 

5.3.1 Symmetric System 93 

5.3.2 Decoupled System 94 

5.3.3 Linear System ' 96 

5.4 Input from The Road Irregularities 97 

5.5 Results and Discussions 98 

5.5.1 Effect of Linear Damping 101 

5.5.2 Effect of Nonlinear Damping 104 

6 2-DOF HEAVE-PITCH MODEL OF A VEHICLE FOR RANDOM 

EXCITATION 112 

6.1 Introduction 112 

6.2 Equations of Motion 113 

6.3 Road Profile Model 114 

6.4 Monte Carlo Simulation 116 

6.5 Results and Discussions 118 

7 CONCLUSIONS 126 

7.1 Conclusions 126 

7.2 Scope of Future Work . 128 

REFERENCES 130 


xi 



APPENDIX 


A Solution of Ordinary Second Order Differential Equations for Three 
Special Non-homogeneous Functions 141 

B Laplace Transform of Equation (5.41) and Its Inverse 145 


Xll 



List of Figures 


2.1 Non-linear base excited system 16 

2.2 The rounded step displacement 19 

2.3 The rounded pulse displacement 20 

2.4 The oscillatory displacement step 22 

2.5 Comparison of analytical and numerical results for the cubic nonlinearity 

in stiffness only. £ = 0.1,7 = 20, e = 0.01, = 1, <5i = 0. O zeroth 

order, * first order, > second order, numerical. 

Curve (a): the rounded step displacement; curve (b): the rounded 
pulse displacement; curve (c): the oscillatory displacement step. ... 33 

2.6 Comparison of analytical and numerical results for the cubic nonlinearity 
in damping only for the rounded step displacement. £ = 0.1, t — 

0.01, £i = 0, 8i = 1 7 = 20 O zeroth order, * first order, > 

second order, numerical 34 

2.7 Comparison of analytical and numerical results for the cubic nonlinearity 
in damping only for the rounded pulse displacement. £ = 0.1, e = 

0.01, 6\ = 0, = 1 7 = 20 O zeroth order, * first order, > 

second order, numerical 35 

2.8 Comparison of analytical and numerical results for the cubic nonlinearity 
in damping only for the oscillatory displacement step. £ = 0.1, e = 

0.01, £i = 0, fix = 1 7 = 20 O zeroth order, * first order, > 

second order, numerical 36 

xiii 



2.9 Comparison of analytical and numerical results for the cubic nonlinearity 
in damping only for the oscillatory displacement step. £ = 0 . 1 , e = 

0.01, = 0, di = 1 7=1 O zeroth order, * first order, t> 

second order, numerical 37 

2.10 The effect of nonlinear cubic damping on the relative displacement for • 

the rounded displacement step. 7 = 50, £ = 0 . 1 , e* = 0 . 

d* = 0, O <5* = 0.01, * d* = 0.05, d* = 0.1 39 

2.11 The effect of nonlinear cubic damping on the relative velocity for the 

rounded displacement step. 7 = 50, £ = 0.1, e* = 0. 

d* = 0, O d* = 0.01, * d* = 0.05, <5* = 0.1 40 

2.12 The effect of nonlinear cubic damping on the velocity for the rounded 

displacement step. 7 = 50, £ = 0.1, e* = 0. d* = 0, O 

8* = 0,01, * 8* = 0.05, 6* = 0.1 41 

2.13 The effect of nonlinear cubic damping on the acceleration for the rounded 

displacement step. 7 = 50, £ = 0.1, e* = 0 . 5* = 0 , O 

5* = 0.01, * d* = 0.05,. 5* = 0.1 42 

2.14 The effect of nonlinear quadratic damping on the acceleration for the 

rounded step. 7 = 50, £ = 0.1, e* = 0. 5* = 0, O 

8* = 0.01, > d* = -0.01 44 

3.1 Non-linear force excited system 46 

3.2 The rounded step force excitation 49 

3.3 The rounded pulse force excitation 50 

3.4 The effect of nonlinear cubic damping on the displacement for the round- 
ed step force excitation. 7 = 50, £ = 0.1. : d* = 

0, .... : d* = 0.05 52 

3.5 The effect of nonlinear cubic damping on the velocity for the rounded 

step force excitation. 7 = 50, £ = 0.1. : d* = 0, .... : d* = 

0.05 53 


xiv 



3.6 The effect of nonlinear cubic damping on the acceleration for the rounded 

step force excitation. 7 = 50, £ = 0.1. : 5* = 0, .... : 5* = 

0.05 54 

3.7 The effect of nonlinear cubic damping on the displacement for the round- 
ed pulse force excitation. 7 = 1, £ = 0.1. : 5* = 

0, .... : 5* = 0.05 55 

3.8 The effect of nonlinear cubic damping on the velocity for the rounded 

pulse force excitation. 7 = 1, £ = 0-1- : 5* = 0, .... : 5* = 0.05 56 

3.9 The effect of nonlinear cubic damping on the acceleration for the rounded 

pulse force excitation. 7 = 1, £ = 0.1. : <5* = 0, .... : 5* = 0.05 57 

4.1 Type 1 Isolator 61 

4.2 Type 2 Isolator 62 

4.3 Type 3 Isolator 64 

4.4 Type 4 Isolator 66 

4.5 SDR and RDR curves for type 1 isolator as a function of 7 for the 

rounded step displacement. I. linear damping £ = 0.1, II. linear + cubic 
damping £ = 0.1,5* = 0.01, III. linear + Coulomb damping £ = 0.1,5/ = 

0.5 and IV. linear + cubic + Coulomb damping £ = 0.1, 5* = 0.01, 5/ = 

0.5 (suffix S for SDR and R for RDR) 68 

4.6 SAR curves in dB for type 1 isolator as a function of 7 for the rounded 
step displacement. I. linear damping £ = 0.1, II. linear + cubic damping 
£ = 0.1, 5* = 0.01, III. linear + Coulomb damping £ = 0.1, 5/ = 0.5 and 

IV. linear + cubic + Coulomb damping £ = 0.1, 5* = 0.01, 5/ = 0.5 . . . 69 

4.7 SDR and RDR curves for type 2 isolator as a function of 7 for the 

rounded step displacement. I. linear damping £ = 0.1, II. linear + cubic 
damping £ = 0.1,5* = 0.01, V. type 2 isolator £ = 0.1,5* = 0.01, n = 1 
(suffix S for SDR and R for RDR) 72 

4.8 SAR curves in dB for type 2 isolator as a function of 7 for the rounded 
step displacement. I. linear damping £ = 0.1, II. linear 4- cubic damping 

£ = 0.1, 5* = 0.01, V. type 2 isolator £ = 0.1, 5* = 0.01, n = 1 73 


xv 



4.9 The effect of stiffness ratio (n) on the SDR in type 2 isolator for the 

rounded step displacement at 7 = 50 (£ = 0.1, 5* = 0.01) 74 

4.10 SDR curves for type 3 isolator as a function of 7 for the rounded step 
displacement. I. linear damping £ — 0.1, II. linear + cubic damping 
C = 0.1,5* = 0.01, VI. type 3 isolator £ = 0.1,5* = 0.01, n = 1 , 5 2 = 

0.01, = 0.5 76 

4.11 SDR curves for type 3 isolator as a function of 7 for the oscillatory 
displacement step. I. linear damping £ = 0.1, II. linear + cubic damping 
£ = 0.1,5* = 0.01, VI. type 3 isolator £ = 0.1,5* = 0.01, n = 1,52 = 

0.01, At = 0.5 . 77 

4.12 SAR curves in dB for type 3 isolator as a function of 7 for the rounded 
step displacement. I. linear damping £ = 0.1, II. linear + cubic damping 
£ = 0.1,5* = 0.01, VI. type 3 isolator £ = 0.1,5* = 0.01, n = 1,52 = 

0.01, ^ = 0.5 78 

4.13 SAR curves in dB for type 3 isolator as a function of 7 for the rounded 
pulse displacement. I. linear damping £ = 0.1, II. linear + cubic damping 
£ = 0.1,5* = 0.01, VI. type 3 isolator £ = 0.1,5* = 0.01, n = 1,52 = 

0.01, = 0.5 79 

4.14 Curve (a): The effect of stiffness ratio (n) on the SDR in type 3 isolator, 
for the rounded step displacement at 7 = 50 (£ = 0.1,5* = 0.01,52 = 

0.01, At = 0.5). Curve (b): The effect of mass ratio (/ 1 ) on the SDR 
in type 3 isolator, for the rounded step displacement at 7 = 50 (£ = 

0.1, 5* = 0.01, 5 2 = 0.01, n = 1.0) 80 

4.15 SDR curves for type 4 isolator as a function of 7 for the rounded step 
displacement. I. linear damping £ = 0.1, II. linear + cubic damping 
£ = 0.1,5* = 0.01, VII. type 4 isolator £ = 0.1,5* = 0.01, n = 1,5 2 = 

0.01, /r = 0.1 82 


xvi 



4.16 SAR curves in dB for type 4 isolator as a function of 7 for the rounded 
step displacement. I. linear damping ( = 0.1, II. linear + cubic damping 
C = 0.1, 6* = 0.01, VII. type 4 isolator ( — 0.1, 6* = 0.01, n = 1,<5 2 = 

0.01, m = 0.1 83 

4.17 The effect of mass ratio (fi) on the SDR in type 4 isolator, for the rounded 

step displacement at 7 = 50 (C = 0.1, 5* = 0.01, 5 2 = 0.01, n = 1.0) ... 84 

5.1 Heave-Pitch Model I 88 

5.2 Heave-Pitch Model II 89 

5.3 Effect of Linear Damping on HDR and PDR for Model I 99 

5.4 Effect of Linear Damping on HAR for Model I 100 

5.5 Effect of Linear Damping on PAR for Model I 102 

5.6 Comparison of Model I and II for HDR and PDR, linear damping case 105 

5.7 Comparison of Model I and II on HAR and PAR, linear damping case 107 

5.8 Effect of Nonlinear Damping on HDR and PDR 108 

5.9 Effect of Nonlinear Damping on HAR and PAR 109 

5.10 Comparison of Model I and II for HDR and PDR, nonlinear damping easel 10 

5.11 Comparison of Model I and II for HDR and PDR, nonlinear damping caselll 

6.1 Comparison of generated autocorrelation function with target autocor- 
relation function. Yi: target; g: generated 119 

6.2 Effect of nonlinear damping on RMS Heave displacement. . . . ■ 122 

6.3 Effect of nonlinear damping on RMS Heave acceleration in dB 123 

6.4 Effect of nonlinear damping on RMS Pitch displacement 124 

6.5 Effect of nonlinear damping on RMS Pitch acceleration in dB 125 


xvii 



List of Symbols 


a 

b 

Co 

Cl 

c 2 
C 3 

ClO 

Cu 

C20 

C21 

c f 

E 

h 

i 

I 

ko 

L 

m 

m 2 

Q- 

R 

Ry 

s 

sgf{x) 


distance from the center of mass to the front axle 
distance from the center of mass to the rear axle 
linear damping coefficient 
nonlinear damping coefficient 

linear damping coefficient of the secondary system 
nonlinear damping coefficient of the secondary system 
linear damping coefficient of the front tyre 
nonlinear damping coefficient of the front tyre 
linear damping coefficient of the rear tyre 
nonlinear damping coefficient of the rear tyre 
Coulomb damping coefficient 
expectation operator 
horizontal lag 

n/=T 

moment of inertia 
linear stiffness constant 
vehicle wheel base 
mass 

secondary mass 
road roughness parameter 
non-dimensional radius of gyration (/><■& 
autocorrelation function of random process y 
horizontal distance variable 
= -1 if x < 0 
= 0 if x = 0 
= +1 if x > 0 




xviii 



Sy 

t 

h 

T 

T v 

T 0 

V 

W(t) 

x 

Vi 

V2 

a 

5 

5* 

Sf 

A 

e* 

fin 

7 

A 

V 

9 

a 

T 

u 

UIq 

n 

C 

kl 

M 

I X I 


V 


Ls-s)) 


*jJi\ Tr\ C 


power spectral density of random process y 
time 

time delay (.*■'*-)) 

non-dimensional time (M*- f ^ « V I 
non-dimensional time period C 1 I 

non-dimensional time delay £ p 3 5b Covv^- 

vehicle forward speed 
white noise process 

vertical displacement of the center of mass 
base excitation due to road irregularities at the front tyre 
base excitation due to road irregularities at the rear tyre 
offset parameter (_ 5 9 <u '-4 9 1/ ^ <m 

Dirac-delta function 

non-dimensional nonlinear damping coefficient ^ ^ 

non-dimensional Coulomb damping coefficient ( 

non-dimensional relative displacement (MJt- Ifi, cr*\ (. % 

non-dimensional nonlinear stiffness coefficient ^ -j 

random phase angle 
shock severity parameter 
wavelength 
mass ratio 
wave number 

angle of rotation about center of mass 
standard deviation of road irregularity 
time lag 

circular frequency 
linear natural frequency 

non-dimensional frequency " y 

non-dimensional linear damping ratio ^ mjl 
largest integer < x 
smallest integer > x 
absolute value of x 


Li &r\ 


( 


UC 


TVf- 


.aM 


*1 


O )})) 

pM'<rr\ C X r 2 '). 


XIX 



Chapter 1 


INTRODUCTION 

1.1 Introduction 

Vibration and shock phenomena are omnipresent. Many occur in nature and many in 
man-made devices and structures. Quite often vibration and shock are not desirable 
and the interest lies in reducing it. Insertion of a vibration isolator between the source 
and receiver of vibration is a very common method of vibration control. The function 
of an isolator is to reduce the magnitude of motion transmitted from the foundation 
to the equipment or to reduce the magnitude of force transmitted from the equipment 
to the foundation. 

Shock is a transient condition. The equilibrium of a system is disrupted either by 
suddenly applied forces or by increment of force or by sudden change in the motion. 
This disruption and its ensuing reaction of the system to restore equilibrium constitutes 
condition of shock. 

Examples of shock are aircraft landing, braking and gust loading, missile launching 
and braking, transportation of fragile equipment, accidental collision of vehicles, drop 
forging machine, vehicle running on a bump, high speed fluid entry, explosions, gunfire, 
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etc. 


The subject of shock is exceedingly complex because the suddenly applied forces 
and sudden motions cause mechanical structures to vibrate in a transient manner, and 
in many natural modes simultaneously. In the application of shock isolators to an ac- 
tual equipment, the location of isolators are determined largely by practical mechanical 
considerations. In general, this results in types of non-symmetry and coupled modes 
not well adopted to analysis by simple means. It is convenient in the design of shock 
isolators to idealize the system as one having symmetry and uncoupled modes of vibra- 
tion. In the simplified analysis the masses are considered to be concentrated and the 
flexible members are massless. The effect of this simplification is to discard the high- 
er modes of vibration. The much less complex analysis of a single degree-of- freedom 
(SDOF) system then becomes applicable. Mounted equipment and isolators assembly 
is represented by a mass-spring-dashpot system. 

There is a considerable justification for this simplification because the cushioning 
afforded by the isolator reduces the suddenness of forces and motions, and thereby 
tends to reduce the response of structures in their higher modes of vibration. 

A different form of idealization is required when the equipment is flexible; e.g., a 
large relatively flexible aircraft subjected to landing shock. Then it is important to 
represent the aircraft as a system with several degrees-of-freedom. 

Isolation of vibration or shock signifies the temporary storage of energy and its 
subsequent release substantially in its entirety but in a different time relation. Isolation 
is thus distinct from the absorption or dissipation of energy. The effectiveness of an 
isolator is sometimes enhanced by limited dissipation of energy but this is a secondary 
consideration in its function. 

A large class of physical systems operate in random vibration environment. These 
include: aerospace systems excited by atmospheric and boundary layer turbulence and 
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jet noise, aircraft and vehicles subjected to track induced vibrations; ground based 
structures excited by earthquakes and wind; and offshore structures excited by wind 
and hydrodynamic wave-induced loads. A probabilistic approach, based on random 
vibration analysis, provides a rational and consistent basis for the analysis and design 
of such systems. 

In the most general case, the excitation can be a combination of transient, periodic 
or random disturbances. It is therefore clear that the optimum design of an isolator 
for a general type of excitation is probably not feasible and in practice an isolator will 
be designed for only one set of environmental conditions. 

The linear shock isolation problems have been widely reported in the literature. 
Nonlinearity is ubiquitous in nature. Linearity is an approximation to reality. In shock 
and vibration systems, isolators such as air springs, elastomeric dampers and wire-rope 
isolators arc inherently nonlinear. Assumption of Hooke’s law for springs and linear 
viscous damping for dampers is done just for mathematical simplicity. Sometimes, 
the amplitudes of steady-state vibration are small enough to justify the assumption of 
linearity. However, the transient displacements may often be sufficiently large when 
the nonlinearity in springs and dampers cannot be ignored. 

The primary objective of this thesis is to analyze the role of nonlinearity in both the 
stiffness and damping on the performance of shock isolation systems. In the present 
work, the isolation system is modelled as SDOF and 2-DOF systems with shock ex- 
citation. At the end, 2-DOF vehicle model has been considered for both shock and 
random base excitations. 



1.2 Review of Previous Work 


The amount of literature available on the general problem of shock and random vi- 
bration isolation is too vast to be considered in detail. In this section, the literature 
pertinent to the present work is discussed. An overview of the research work in dif- 
ferent areas, connected with the present problem, is categorized under the following 
headings: 

(i) Shock Isolators, 

(ii) Random Vibration Isolators 

1.2.1 Shock Isolators 

Linear shock isolation problems are discussed in several books [1-4]. Snowdon [4] 
presented the response of a nonlinear shock isolator modelled as a nonlinear elastic 
(tangent and inverse tangent elasticity) spring parallel to a viscous damper. He con- 
cluded that a soft spring (inverse tangent elasticity) performs better than a hard spring. 
In reference [5], Snowdon compared the performance of a dual-phase damper mounting 
system with that of a linear, simple mounting system. Guntur and Sankar [6] reported 
the performances of different kinds of dual-phase damper mounting system and com- 
pared these with that of a linear, simple mounting system. Hundal [7] reviewed the 
literature on pneumatic shock absorbers and isolators. Hundal [8] has also compared 
the performances of shock isolators with linear and quadratic damping with a base 
input in the form of an acceleration pulse of rectangular shape. 

By using the concept of a variable friction force, Mercer and Rees [9] proposed a new 
form of shock isolator which is adaptive in its action but is still composed of entirely 
passive elements. Eshleman and Rao [10] have investigated shock isolation characteris- 
tics of six types of elements: (a) helical coil spring, (b) ring spring, (c) friction snubber, 
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(d) liquid spring, (e) pneumatic spring, and (f) solid rubber elastomer. They concluded 
that the friction snubber-helical spring combination yielded the best isolation of all the 
elements. Hundal [11] describes the analysis of response of a pneumatic shock isolator 
to base acceleration of rectangular and half-sine shape pulses. The isolator consists of 
a pneumatic dampei in parallel with a linear spring. Hundal [12] has also presented 
the analysis and lesponse of a symmetric pneumatic absorber with orifice damping. 
The modelling and response of a pneumatic shock isolator are presented in reference 
[13]. The isolator consisted of a symmetric pneumatic spring with self-damping and 
was subjected to a velocity pulse. In reference [14], an unsymmetric pneumatic shock 
isolator is considered. An investigation on a linear shock isolator subjected to a velocity 
pulse is included in reference [15]. 

The concept of an ‘on-off’ damper employing the feedback signals from directly mea- 
surable variables is proposed in reference [16]. A control scheme utilizing the directly 
measurable relative position and relative velocity signals to produce the command sig- 
nal is configured. For larger values of shock severity, the ‘on-off’ dampers perform 
extremely superior to any passive damper. 

The shock isolation characteristics of semi-active isolators based on two different 
control schemes are presented in reference [17]. In the first type, the control depends 
on the absolute and relative velocities of the sprung mass while in the second type it 
depends on the relative displacement and the relative velocity. Both employ continuous 
control of damper forces as opposed to ‘on-off’ control. The damper forces are generated 
by modulating its orifice areas for fluid flow. The constraint on damper force is that 
the energy associated with it is always dissipated. It is seen that both the semi- 
active control schemes can simultaneously reduce the acceleration and the relative 
displacement, something passive linear isolators cannot do. 

Various analytical methods are available in the literature [18-32] for transient analy- 
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sis of nonlinear systems, namely, the linearization method [19], the ultraspherical poly- 
nomial method [20-22] and Lighthill’s extension of Poincare’s perturbation method [23]. 
A simple-straight forward perturbation method is developed in [32] by using Laplace 
transforms. 

1.2.2 Random Vibration Isolators 

Numerous papers have so far been published on random vibration isolators especially 
with reference to vehicle vibration control. In reference [33] various vehicle models 
have been presented. Hrovat [34] studied the preview control applied to a linear, two- 
dimensional (2-D), four-degree- of-freedom, half-car model. The model includes vehicle 
sprung mass pitch and heave dynamics and front and rear unsprung masses. Linear 
Quadratic (LQ) control formulation has established that the preview based on front 
wheel input can contribute up to 70% reduction of vehicle r.m.s acceleration, for a given 
suspension travel, and up to 50% reduction in tire deflection for a given acceleration 
level. In reference [35], a fully active fail safe suspension is reported for the 4-DOF 
model using the optimal regulator theory. Root mean square (r.m.s) bounce and pitch 
absolute acceleration transmissibility of the sprung mass and of maximum suspension 
and tire deflections are presented for semi-active suspensions. Both passive and active 
schemes have been considered. In reference [36], the response of a vehicle moving over 
a non-homogeneously rough flexible track supported over a complaint foundation has 
been obtained for variable and constant velocity runs. The characteristics of the vehicle 
and track mode responses are expressed in closed form for both heave-pitch and heave- 
roll models. A quarter car suspension system model, containing a controllable damper 
with a limited range of coefficient values, a limited adjustment system bandwidth and a 
realistic elastic mounting to the car body has been optimized for operation on a random 
road [37]. The optimization has been carried out through a novel numerical method. In 
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reference [38], an approximate state-space method for obtaining the time varying mean 
and covariance of the response of nonlinear systems excited by non-stationary random 
process is presented. Redfield and Karnopp [39] analyzed a two-degree-of-freedom 
model for its performance sensitivity as a function of the system’s feedback gains. A 
two-degree-of-freedom (1/4 car model) has also been used to evaluate alternative linear 
control laws [40]. Hac [41] considered decentralized control of active vehicle suspensions 
with preview of road irregularities using a two-degree-of-freedom vehicle model. He has 
also studied the problem of active suspension control of a two-degree-of- freedom vehicle 
traveling on a randomly profiled road [42]. 

Kirk [43] has investigated the influence of nonlinear spring stiffness characteristics 
on the effectiveness of vibration isolators with linear damping, subjected to stationary, 
random white noise ground acceleration. The probability density function of spring 
displacement is determined analytically by means of the Fokker-Planck equation and 
both r.m.s and mean peak values of spring displacement and mass acceleration are pre- 
sented for three types of spring nonlinearity: (a) cubic hard spring, (b) cubic soft spring 
and (c) tangent spring. A method for modelling vehicles as a system of interconnected 
rigid and flexible bodies is presented in reference [44]. Finite element and component 
mode substitution techniques are employed to characterize elastic motions of flexible 
bodies in terms of minimal set of modal co-ordinates. Semi-active control techniques 
are then used for suppressing vibrations at some desired points, in the chassis, located 
away from suspension attachments. In references [45-46], a preview control problem 
of a semi-active suspension was formulated and a truly optimal solution was obtained 
using the calculus of variations. 

Crolla and Abdel-Hady [47] have studied the effects of the following issues in de- 
riving control laws based on a mathematical model of an actively suspended vehicle: 

(a) representation of the ground surface as integrated or filtered white noise, 
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(b) cross-correlation between left and right track inputs, 

(c) wheelbase time delay between front and rear inputs. 

In reference [48], the basic functions which suspensions must perform are described 
and categories of controllable suspensions based on the required amount of control 
power are defined. The kinds of devices which can be used in controllable suspensions 
are discussed as well as several basic philosophies. Finally, fundamental theoretical and 
practical limitations in suspension performance are indicated. Lin and Willumeit [49] 
have simulated the dynamics of a vibrator-controlled adaptive damper. Its potential in 
vehicle isolation improvement is demonstrated by a quarter car model. The problem of 
deriving control laws which minimize specified performance indices for a vehicle moving 
on a rough surface with preview of the surface elevation is considered in reference [50]. 
Bhave [51] has presented an analysis of a vehicle mounted on air springs at the front 
and rear which are connected through a capillary tube. 

Fodor and Redfield [52] have introduced the concept of regenerative damping in vi- 
bration control, the storage of energy normally dissipated by a passive viscous damper. 
A device, called the variable linear transmission, is proposed to accomplish this task. 
In reference [53], a continuous time control strategy for an active suspension with pre- 
view, based on optimal control theory is presented. It has been shown in reference [54] 
that an electromagnetic rotating damper can be used as a controlled force generator 
for vibration isolation in vehicles. Moran and Nagai [55] have shown that how neural 
networks can be used for the forward and inverse identification as well as the optimal 
control of nonlinear active suspensions. Well trained neuro-vehicle models can identify 
the dynamics of actual vehicle suspensions and predict their future behaviour. 

Narayanan and Raju [56] have investigated an active suspension system to control 
the non-stationary response of a single-degree-of-freedom (SDOF) vehicle model with 
variable velocity traverse over a rough road. They considered the design of an active 
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optimal suspension to control the non-stationary response of a 2-DOF vehicle model 
traversing a rough road [57, 58], They have also studied the active control of non- 
stationary response of a two-degree-of-freedom vehicle model with nonlinear passive 
suspension elements [59]. The method of equivalent linearization is used to derive 
an equivalent linear model and optimal control laws are obtained by using stochastic 
optimal control theory based on full state information. Velocity squared quadratic 
damping and hysteretic type of stiffness nonlinearities are considered. In reference 
[60], active control of the time varying response to a stationary random excitation 
of a two degree-of-freedom (2-DOF) vehicle model with nonlinear passive suspension 
elements have been considered. 

Tamboli and Joshi [61] have obtained the Root Mean Square Acceleration Response 
(RMSAR) of a vehicle dynamic system subjected to actual random road excitations, so 
as to account for the effect of the actual power spectral density (psd) of road excitation 
and the frequent changes in vehicle velocity. A statistical dimensionless ratio, a ratio 
of excitations at rear and front axles, has been formulated to study the effect of time 
lag on the response of the vehicle due to wheelbase. The effect of vehicle velocity and 
wheelbase on the RMSAR has been studied and the RMSAR is found proportional to 
the vehicle velocity and inversely proportional to the wheelbase. 

Ariaratnam and Wei-Chau Xie [62] have presented various methods to investigate 
the almost-sure asymptotic stability of nonlinear coupled oscillators under external 
random excitation. In reference [63], an asymptotic expansion for the maximal Lya- 
punov exponent, the exponential growth rate of solutions to a linear stochastic system, 
and the moment Lyapunov exponent, the asymptotic growth rate of the moments of 
the response, have been obtained for systems driven by a small intensity real noise 
process. The results obtained are applied to a thin rectangular beam under the action 
of a stochastic follower force and a model of a vehicle traveling over a rough road. 
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Blevins [64] has reported the analysis to determine the properties of a random 
process consisting of the sum of a series of sine waves with deterministic amplitudes 
and independent random phase angles. The probability density of the series, its peaks 
and envelope have been found for an arbitrary number of sine waves in the series. 

Shinozuka [65] has presented the various efficient methods for digital simulation 
of a general homogeneous process (multidimensional or multivariate or multivariate- 
multidimensional) as a series of cosine functions with weighted amplitudes, almost 
evenly spaced frequencies and random phase angles. The approach is also extended to 
the simulation of a general non-homogeneous oscillators process characterized by an 
evolutionary power spectrum. Simulation of one-dimensional uni-variate, stationary, 
Gaussian stochastic processes using the spectral representation method was presented 
in reference [66]. In reference [67], a stochastic model for general periodic excitations 
with random disturbance which is constructed by introducing random amplitude and 
phase disturbances to individual term in the Fourier series of the corresponding de- 
terministic periodic function has been presented. Johnson et. al.[68] compared the 
performance of Monte Carlo simulation on various platforms, in the context of a four- 
dimensional linear oscillator and a Duffing oscillator subjected to a band-limited white 
noise. 

Zeldin and Spanos [69] presented the representation and simulation of random fields 
using wavelet bases. Examples of simulating random fields encountered in engineering 
applications are discussed. A Haar wavelet based model for probabilistic behaviour of 
stochastic dynamic systems was presented [70]. A Hamilton’s law of varying action 
based approach was used to obtain a deterministic model. A single degree-of- freedom 
spring-mass-damper system subjected to a random forcing function was considered 
to show the feasibility and the effectiveness of the formulation. In reference [71] a 
new hierarchical method for the Monte Carlo simulation of random fields called the 
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Fourier-wavelet method is developed and applied to isotropic Gaussian random fields 
with power law spectral density functions. This technique is based upon the orthog- 
onal decomposition of the Fourier stochastic integral representation of the field using 
wavelets. 


1.3 Objectives and Organization of the Present Work 

The primary objective of the present thesis is to study the effect of nonlinearity in the 
isolator elements subjected to shock excitation. Both base and force shock excitations 
i have been studied. A simple analytical tool has been developed using perturbation 

method in conjunction with Laplace transform to determine the transient response of 
an SDOF system in the presence of nonlinear, dissipative shock isolators. To improve 
the performance of a shock isolator with base excitation, in the presence of nonlinear 
cubic damping, four methods namely, Coulomb damped isolator, three element isolator, 
vibration absorber and two stage isolator have been compared. Next, a 2-DOF (Heave 
and Pitch) vehicle model moving with constant speed, with a time delay in front and 
rear wheel inputs, has been studied for shock and random excitations. The organization 
of the thesis is detailed below. 

In chapter 2, a simple technique combining the straightforward perturbation method 
with Laplace transform has been developed to determine the transient response of a 
single degree-of-freedom system in the presence of nonlinear, dissipative shock isolators. 
Analytical results are compared with those obtained by numerical integration using 
the classical Runge-Kutta method. Three types of input base excitations, namely, the 
rounded step, the rounded pulse and the oscillatory step are considered. The effects 
of nonlinear damping on the response are discussed in detail. Both the positive and 
negative coefficients of the nonlinear damping term have been considered. It has been 
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shown that a critical value of the positive coefficient maximizes the peak values of 
relative and absolute displacements. This is true for any power-law damping force 
with an index greater than 1. On the other hand, the overall performance of a shock 
isolator improves if the nonlinear damping term is symmetric and quadratic with a 
negative coefficient as exhibited by elastomers up to a certain value of the amplitude of 
deformation. It has been established that the nonlinearity in damping predominantly 
affects the shock isolator response, whereas the effect of nonlinearity in stiffness is 
insignificant. So, in the subsequent chapters only the damping (and not the stiffness) 
is taken to be nonlinear. 

The performance of the isolator system considered in chapter 2 is re-examined 
under force shock excitations. These results are included in chapter 3. Unlike with 
base excitation, the presence of nonlinearity in the damping force is seen to marginally 
improve the performance of a linear isolator if the system is excited by a shock force 
on the mass. 

In the presence of inherent nonlinear cubic damping in a nonlinear shock isolator 
with base excitation, four different modifications of the isolator are considered to im- 
prove its performance. These modifications are (i) an isolator with a Coulomb damper, 
(ii) a three element isolator with the nonlinear damper mounted elastically, (iii) an iso- 
lator with vibration absorber and (iv) a two stage isolator. Overall, it is seen that the 
three element and two stage isolators are preferable in the presence of nonlinear cubic 
damping. These results are presented in chapter 4. 

In chapter 5, a 2-DOF model for a vehicle undergoing heave and pitch motions 
is investigated for shock excitations. It is assumed that the vehicle is moving with a 
constant speed so that the rear wheels are subjected to the same input as the front 
wheels but with a time delay. Both coupled and decoupled symmetric systems are 
considered. It is shown that the decoupled symmetric system gives optimum isolation 
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performance. Like in SDOF systems, nonlinear damping is detrimental to 2-DOF 
systems for shock excitation with a two element isolator. A three element isolator 
nullifies the adverse effect of nonlinear damping in 2-DOF system as well. 

The vehicle model considered in chapter 5 is investigated again with random ex- 
citation. This study is presented in chapter 6, Here, the front and rear suspensions 
are subjected to the same filtered white noise but with some time delay governed by 
the speed of the vehicle. Shinozuka's [66] spectral method is used for generation of the 

random process. Monte Carlo simulation has been employed to find the response of 
the vehicle. 

In chapter 7, overall conclusions and scope of future work have been presented. 
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Chapter 2 


IESPONSE OF NONLINEAR 
DISSIPATIVE SHOCK 
[SOLATORS TO 

DISPLACEMENT EXCITATION 

2.1 Introduction 

Nonlinearity is ubiquitous in nature. Linearity is an approximation to reality. In shock 
and vibration systems, isolators such as air springs, elastomeric dampers, and wire-rope 
isolators are inherently nonlinear. Assumption of Hookes’ law for springs and linear 
viscous damping for dampers is done just for mathematical simplicity. Sometimes the 
amplitudes of steady-state vibration are small enough to justify the assumption of 
linearity. However, the transient displacements may often be sufficiently large when 
the nonlinearity in springs and dampers cannot be ignored. 

Ravindra and Mallik [72-74] have studied the effect of nonlinear damping on the 
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performance of vibration isolators under harmonic loading. They observed bifurcations, 
chaos and strange attractors due to the presence of nonlinearity in springs and dampers 
in vibiation isolators. They concluded that a strictly dissipative nonlinear damping 
may be used as a passive conti ol strategy to suppress various instabilities occurring 
in nonlinear vibration isolation systems. The major objective of this chapter is to 
ascertain the effect of nonlinear damping on the response of shock isolators. 

Vaiious analytical methods are available in the literature [18-32] for transient analy- 
sis of nonlinear systems, namely, Linearization method [19], Ultraspherical polynomial 
method [20-22], and Lighthill’s extension of Poincare’s perturbation method [23], In 
this chapter, a simple straightforward perturbation method along with Laplace trans- 
form is used. This method is applicable for any order of nonlinearity, both in the 
restoring and damping forces, expressed in the form of polynomials. A cubic non- 
linearity in the restoring force is assumed, whereas various forms of nonlinearities so 
far as the damping force is concerned have been included. For example, a cubic type 
nonlinearity over and above the common linear viscous damping, as exhibited by fluid- 
dampers at high velocities has been given special attention. Numerical results are 
included for a typical elastomeric damper which can be modelled up to a certain am- 
plitude of deformation, by a combination of linear viscous damping and a dissipative 
quadratic damping with a negative coefficient[75]. 

Three types of base excitations, namely, the rounded step, the rounded pulse and 
the oscillatory step given in reference [5] are considered. Analytical results obtained by 
the present method are compared with those obtained by direct numerical integration. 
For the quadratic damping case, results are obtained by numerical integration alone 
since the symmetric dissipative quadratic damping can not be expressed in polynomial 
form. The effects of nonlinear damping on the response of various models of shock 
isolators are discussed in detail. 
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2.2 Equations of Motion 

A single degree-of-freedom (SDOF) shock isolator system is shown in Figure 2.1, where 
the base is subjected to a shock displacement. The spring and the damper elements 
of the isolator are taken to be nonlinear with a cubic nonlinearity superimposed on a 
linear term. The equation of motion for the mass m, to be isolated, is 

mx 2 + c 0 (±2 - ±0 + ci (± 2 - xxf + k 0 {x 2 - Xl ) + ki(: r 2 - x x ) 3 = 0 (2.1) 

where X\ and x 2 are the absolute displacements of the base and the mass m, respectively, 
and the dots denote derivatives with respect to time t, fco and k\ are linear and nonlinear 
stiffness constants, and Co and c x are linear and nonlinear damping coefficients. The 
initial conditions on x 2 are taken to be zero, i.e., x 2 = x 2 = 0 at t = 0. This equation 
is written in the non-dimensional form as 

A" + 2CA' + <$*( A') 3 + A + e*A 3 = f(T) (2.2) 

where, A = (x 2 — X\)/x\ max is the relative displacement, between the mass m and the 
base, non-dimensionalised with respect to the maximum base displacement Xi max , T is 
the non-dimensionalised time T = cj 0 t, where w 0 = is the linear natural frequen- 
cy, the primes denote derivatives with respect to T, and the other non-dimensional 
parameters are 



/(T) = - Ah. (2.3) 

*£l max 

2.3 Types of Input Displacements of The Base 

Three different input displacements [4, 5] considered in this study are defined in the 
following section. These inputs are chosen to be physically realistic in that they describe 
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the translation of the foundation through a finite distance in a finite time with finite 
acceleration and deceleration. These inputs are defined by the following equations: 

Case (a) Rounded displacement step: 

This type of base excitation is shown in Figure 2.2 and is described as 
Xi(t) = 0 when t < 0, 

xi(t) = X\ max 1 - (1 + 7w 0 t) e (_7 “ o4) ] for t>0 (2.4) 

with 7 as the severity parameter. The severity parameter (7) describes the finite rise 
time, r, of the step in terms of the half-period of natural vibration T v j 2 of the isolator 
system; thus 

7 = T p /(2r) = 7r /(lvqt). (2.5) 

The rise time r of the rounded displacement step is defined as the time required for 
the displacement to reach 82% of its final value. 

Case (b) Unidirectional rounded displacement pulse: 

This type of pulse excitation is shown in Figure 2.3 and is described as 

Xi(t) = 0 when t < 0, 

Xi(t) = x lmax (e 2 /4) (7w 0 i) 2 e (_7a,ot) for t > 0. (2.6) 

Again, the severity parameter (7) is defined by equation (2.5). Here, the duration r 
of the displacement pulse is defined as the length of an equivalent rectangular pulse 
that has the same area as that of the rounded pulse, but with a height 17.6% 

higher than £ lmai . 
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Case (c) Oscillatory displacement step: 

This type of base excitation is shown in Figure 2.4 and is described as 
Xi(t) = 0 when t < 0 , 

= X imax (0.68684) j^l — {cos( 7 o;of) + 0.25 sin(7ajot)}e^""°' 257 ‘ Joi )j (2.7) 

for t > 0 . 

Once again, the severity parameter ( 7 ) is defined by equation (2.5), but here, the 
rise time r of the oscillatory displacement step is exactly the time required for the 
displacement to reach its maximum value x lmax . 

The method of selecting the input displacement by equation (2.7) is explained in 
reference [5]. For the sake of completeness, the method is presented below. This input 
represents the character of an abrupt transient disturbance after it had passed through 
and had been ‘filtered’ by the foundation structure that supports the mounting system 
under consideration. 

The base of an SDOF system is subjected to a step displacement x(t ) = xq , the 
resultant displacement x\ of the system mass is given in reference [4] by -the equation 

x x = £ 0 [ 1 . + e^ FWFt) {(( F UF/ud)sin(ijo d t) - cos(u d t)}] ( 2 . 8 ) 

where ujp is the natural frequency of the filtered system, £f is its damping ratio, and 
u>d = u)F\jl — Cf- 

This input displacement x\ in equation ( 2 . 8 ) is not physically realistic because its 
velocity is discontinuous and its acceleration is infinite at the time origin. However, if 
the sine term in equation ( 2 . 8 ) is reversed in sign, then velocity is zero and acceleration 
remains finite at the time origin. Further, for all but large values of £f> the equation 
may be simplified as follows: 

x'i = Xq[1 — e^ FWF ^{(CF)sin(wp£) + cos(wpf)}]. (2-9) 


21 





This displacement takes a maximum value 


Xlmax — 3?o[l T 6^ ^ F7r ^] 


( 2 . 10 ) 


when 


u F t = 7 r (2.11) 

this equation directly specifies the rise time of the step. Consequently, from equation 
(2.5) it is possible to write that 


7 = < 2 - 12 ) 

and, from equation (2.9) one can write 


^ ^ ^ 0 < HCFsin( 7 u;ot) + cos( 7 o; 0 f )}]/[! + ^ F7r ^]. (2.13) 


The oscillatory displacement step of equation (2.7) is obtained by substituting Q = 0.25 
in equation (2.13). 

The maximum values attained by the velocity and acceleration of the above three 
input displacements of equations (2.4), (2.6) and (2.7) are listed here, because these 
values will be referred in subsequent chapters when the various performance indices of 
shock isolators are evaluated. 

For rounded step 


Xlrnax — 0.367887^0^1 max 3 

__ 2 2 

X\max — T ^Q^lmax 1 

for rounded pulse 

Xi max = 0. 8518870)0^1 max: 

m ax = 3.694537 2 cj 0 Xi maa: 


(2.14) 

(2.15) 


(2.16) 

(2.17) 
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and for oscillatory step 


Xi max = 0.508257WOS1 max 5 (2.18) 

Ximax = 0 . 729777 2 Wq £ lmax • (2.19) 

2.4 Closed Form Solutions 

It is desirable to have an analytical solution of equation (2.2). This would facilitate 
calculating transients at any given instant directly from the closed form solution rather 
than integrating numerically equation (2.2) from T = 0. To this end, perturbation 
method [23-32] has been used to separate the terms of the same order and then the 
solutions of differential equations are obtained through Laplace transformation. 

The parameters 5* and e* in equation (2.2) are taken to be small and of the same 
order, so that equation (2.2) can be rewritten as 

A" + 2(A' + A + e[5i(A') 3 + e x A 3 ] = f(T) (2.20) 


where e<l. 

The form of equation (2.20) is useful in that, by putting e\ and 5i equal to zero or 
one, the effect of nonlinearity, in only stiffness or in only damping or in both stiffness 
and damping, can be studied. The solution of equation (2.20) is perturbed in small 
parameter e(<C 1) as 

00 

A(T) = ]TVAi(T). (2.21) 

i = o 

Substitution of equation (2.21) in equation (2.20) and collection of like power terms of 
e yields 

A" + 2CA' 0 + A 0 = f{T) for * = 0 (2.22) 
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and 


A " + 2£A' + Aj — — Cj 

l — l j 


E E 


j=Q n=0 


j = 0 n=0 


(2.23) 


for i > 1, 


with Aj = A( = 0 at T = 0 V i > 0. 


Convergence 

The series in equation (2.21) is a Cauchy product of £“o ei and ££o Aj(T). 
For e < 1, £T 0 eI is a convergent geometric series and £“ 0 A j(T) is bounded 
for ( > 0, 7 > 0 since the solution of equations (2.22) and (2.23) contain exponential 
terms with negative indices, which are bounded. Therefore, the series £g 0 e l ’A j(T) is 
convergent [76]. It is to be noted here that the damping term is retained in equation 
(2.22) in order to avoid secular terms and thereby ensuring bounded solutions [25]. 

The focus here is on transients alone and that too over a short time duration. It was 
observed that the three term solution of equation (2.20) (i.e., up to i = 2, in equation 
(2.21)) give good convergence. Hence the solution of equation (2.20) are given below 
up to second order (i.e., i = 2 in equation (2.21)) only. 

2.4.1 Zeroth Order Solution 

With the initial conditions A = A' = 0 at T = 0, the solutions of equation (2.22), for 
the three input functions in equations (2. 4), (2. 6) and (2.7), are obtained using Laplace 
transformation and are listed below: 
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Case (a) Rounded displacement step: 


Substitution of equation (2. 4) in equation (2.3) yields 

f(T) = - 7 2 (1 - tT) e ( ~ 7r) for T > 0. 


(2.24) 


The equation (2.24) may be rewritten as 

N 

I 

\m — 0 


f(T) = e aT ( £ e ( ~ 7T) for T>0 


(2.25) 


where, a — — 7 , N = l, A 0 = — 7 s , A\ - 7 3 . 

The solution of equation (2.22) for f(T) given in equation (2.25) is obtained by using 
the formulae listed in Case (i) of Appendix A. Therefore, the zeroth order solution may 
be written as 


where, 


MT) 

— Aoi(T) + Ao2(T) 

(2.26) 

Aoi(T) 

= e (_cr) [P! cos (pT) + Qi sin (J3T)]/P\ 

(2.27) 

A 02 (T) 

= e^ aT) [B 1 + B 2 T}-, 

(2.28) 

P 

= vT^T); 

(2.29) 

Qi + iPx 

A\ + Aq (c^ — a) 

(a — a) 2 

(2.30) 

a 

— + iP\ 

(2.31) 

i 

= 

(2.32) 

Bi 

= [a, - 2 (c - 7 )b 2 ) /r ; 

(2.33) 

b 2 

= Ai/r ; 

(2.34) 

r 

= 7 2 - 2qC + 1- 

(2.35) 
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Case (b) Unidirectional rounded displacement pulse: 


By substituting equation(2.6) in equation (2.3), f(T) is obtained as 

/(T) = -(e 2 7 2 / 4 ) (2 - 4 7 T 1 + 7 2 T 2 ) for T > 0. (2.36) 

The equation (2.36) may be rewritten succinctly as 

/CO = e“ r f E e ( -^ for T > 0 (2.37) 

\m= 0 / 

where, a = - 7 , iV = 2, ^4 0 = -e 2 7 2 / 2 ,A! = e 2 7 3 ,A 2 = -e 2 7 4 / 4 . 

The solution of equation (2.22) for /(T) given in equation (2.37) is obtained by 
using again the formulae listed in Case (i) of Appendix A. Therefore, the zeroth order 
solution may be written as 


A 0 (T) = 

A 01 (T) + A 02 (T) 

(2.38) 

Aoi(T) = 

c 0 s(^T) + Qi sin(/?T)]//3; 

(2.39) 

A 02 (T) = 

e (aT) [B x + B 2 T + B 3 T 2 / 2] ; 

(2.40) 

0 = 

7(i - < 2 ); 

(2.41) 

h + ^1 = 

2y4.2 4i(cb — g) 4" 4.0 (a — 

(2.42) 

(a - a) 3 

= 

~C + ifii 

(2.43) 

£ = 

vTTi 

(2.44) 

Bi = 

[ 4 , - 2 (c - 7 )b 2 - b 3 ) /r, 

(2.45) 

b 2 = 

CO 

I 

CQ 

1 

(2.46) 

Bz = 

A 2 /r ; 

(2.47) 

r = 

7 2 - 27 C 7- 1 . 

(2.48) 
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Case (c) Oscillatory displacement step: 


Substitution of equation(2.7) in equation (2.3) yields 


f(T) = -(17/16)(0.68684)7 2 [cos( 7T) - 0.25 sin(7T)]e ( “°- 257r) (2.49) 


for T > 0. 

The solution of equation (2.22) for f(T) given in equation (2.49) is obtained by using 
the formulae listed in Case (ii) of Appendix A. Therefore, the zeroth order solution 
may be written as 


where, 

and 


AoCO 

= 

Aoi(T)+A 02 (T) 

(2.50) 

Aoi(T) 

— 

e (-cr)[p 2 cos(/?T) + q 2 S in (/3T))//3; 

(2.51) 

A«CT) 

= 

e^lRi cos(loiT ) + S\ sin(oo 1 T')]; 

(2.52) 

e 

= 

\/(i - < 2 ); 

(2.53) 

a 

= 

—( + i/3; 

(2.54) 

i 


7(-i); 

(2.55) 

Q 2 + ^2 

= 

Ci(jCt — Ci) “b D 1 UJ 1 
(o — a) 2 + w 2 

(2.56) 

Qfi 

= 

a + iuii ; 

(2.57) 

Si 4" iR\ 


D\ 4- iC\ 

(2.58) 


CK 2 + %C a l "b 1 

Ci 

= 

— (17/16) (0.68684) 7 2 ; 

(2.59) 

D 1 

rr 

(17/16)(0.68684)7 2 /4; 

(2.60) 

a 

= 

—7/4; 

(2.61) 

uq 

= 

- 7 / 4 . 

(2.62) 
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2.4.2 First Order Solution 

For i = 1, equation (2.23) yields 

A" + 2CA' x + Ax = - [e^Ao) 3 + <5 i(Aq) 3 ] . (2.63) 

Case(a) Rounded displacement step: 

From equation (2.26), the non-homogeneous term Aq in equation (2.63) may be 
expressed as 

Ao (T) = Aq X (T) + A 3 2 (F) + 3Ag 1 (r)A 02 (T) + 3A 0 i(T)Ag 2 (r). (2.64) 

From equation (2.27), Ag,(T) in equation (2.64) may be written concisely as 
&oi( r ) = e’ lT [y){C„cos 

_m=l 

where, a x = -3(; wi = /3; w 2 = 3/3; N = 2; Ci = 3Pi(Px 2 + <Qf)/4/5 3 ; 

A, = 3Qi(Pi 2 + Q?)/4/3 3 ; C' 2 = P 1 (P 1 2 -3Q2)/4 i d 3 ; and D 2 = -Qi(Q? - 3P 2 )/4/3 3 . 

From equation (2.28), the second term in the right-hand-side of equation (2.64) 
may be written in simplified form as 

A 3 2 (T) = e° 2T [f)24 J7i r m ] - (2.66) 

.772=1 

where, a 2 = —37; N = 3] Ao = Bf\ A x = 3 B\B 2 \ A 2 = 3B X B 2 and A 3 = B 2 . 

From equations (2.27) and (2.28), the third and the fourth terms in the right-hand- 
side of equation (2.64) are expressed in succinct form as 

- l v 

3Ao X Ao 2 = e azT [C + E cos(FlT) + Fsin(ftT)] £ A m T m _ (2.67) 

|_771=0 

where, a 3 = —(2( + 'y); C = (P 2 4- Qi)/2; E = (P x - Q\)/ 2; P = PiQi! 

0 = 2/3; N= 1; A 0 = 3CB l /P 2 and A r = 3CB 2 /p 2 


and 

" N 



3A 01 Ag 2 = e a * T [Fcos(QT) + Fsin(QT)] ]T A m T m 

(2.68) 


Lm=0 



(w m T) + Dm sin(6 o m T) } (2.65) 
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where, = — (C+27); E — P\\ F = Q\\£l — N — 2 \ A$ — B\\ A\ — 2B1B2 and A2 


Bl 

The sum of the all four terms in the right-hand-side of equation (2.64) can be 
expressed as the sum of the following three generalized functions with appropriate 
values of coefficients and indices 


,oT 


N 

.771=0 


r n 




Yj {Crn COS (w m T) + D m sin(w m T)} 

Lm=0 


and 


e aT [E cos(flT) + F sin(fiT)] 


AT 


E 


Lm=0 


(2.69) 

(2.70) 

(2.71) 


Similarly, the second non-homogeneous term in equation (2. 63) can also be expressed 
as the summation of three generalized functions given by equations (2.69), (2.70) and 


(2.71) . 

The solution of equation (2.63) is obtained by using the formulae listed in Appendix 
A cases (i), (ii) and (iii), corresponding to the terms in equations (2.69), (2.70) and 

(2.71) respectively. 

Case(b) Unidirectional rounded displacement pulse: 

- In this case also, the non-homogeneous term in equation (2.63) can be expressed as 
the summation of the three generalized functions given by equations (2.69), (2.70) and 

(2.71) . Again the solution of equation (2.63) is obtained by using the formulae listed 
in Appendix A. 

Case(c) Oscillatory displacement step: 


In this case, the non-homogeneous term in equation (2.63) can be expressed as 
the generalized function given by equation (2.70). The solution of equation (2.63) is 
obtained by using the formulae listed in case(ii) of Appendix A. 
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2.4.3 Second Order Solution 

For i = 2, equation (2.23) yields 


A" + 2CA'+A 2 = -3[e 1 (A 0 ) 2 A 1 +<5 I (A') 2 A , 1 ]. (2.72) 

Cases (a) and (b): In these two cases, the non-homogeneous term in the equa- 

tion (2.72) can be expressed as the summation of three generalized functions given in 
equations (2.69), (2.70) and (2.71). Once again, the solutions of equation (2.72) are 
obtained by using the formulae given in Appendix A. 

Case (c): In this case, the non-homogeneous term in the equation (2.72) again turns 
out to be only the second generalized function given by equation (2.70). The solution 
of equation (2.72) is obtained by using the formulae listed in case(ii) of Appendix A. 


2.5 Results and Discussions 


The results are presented in two parts: the first one to validate the accuracy of the 


closed form solutions and the second part to highlight the effects of nonlinearity in 
isolator elements. The linear damping ratio ( is taken to be 0.1 in all the results. 

2.5.1 Validation of Closed Form Solutions 


The closed form solutions of equation (2.2), obtained through Laplace transformation of 
equations (2.22) and (2.23), were checked against the numerically integrated solutions 
of equation (2.20). Numerical integration was carried out by using classical fourth 
order Runge-Kutta method. F»e , o u ^ ^ 

* .-o 

with mr^ ■ . , t i 

ICcl- wcrwl'iVe^ -tCe. . ; 
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First the nonlinearity in stiffness alone is considered by taking e = 0.01, e x = 
1 and = 0 in equation (2.20). Figure 2.5 shows the comparative results for the 
three types of excitation, where curves (a), (b) and (c) correspond to three different 
excitations, described by equations (2.4), (2.6) and (2.7) respectively. It was observed 
that for the lower values of 7 < 20 the effect of nonlinearity is not predominant both 
in the nonlinear damping and nonlinear stiffness cases, but for higher values of 7 > 
20 nonlinear damping effect is very high. So, the severity parameter of the shock 
displacement, 7 , is taken to be equal to 20 for all these results. It is clear from 
these curves that even the zeroth order solutions given by equations (2.26), (2.38) and 
(2.50) are very close to the exact numerical solutions. This indicates that the effect of 
nonlinearity in stiffness of the isolator is negligible. The same behaviour was observed 
for a large range of the severity parameter 7 up to 100 . 

The closed form solutions for isolator with nonlinearity in damping alone are pre- 
sented for three types of base excitation in Figures 2.6, 2.7 and 2.8, with e = 0.01, e x = 
0, d x = 1 and 7 = 20. It is evident from these figures that the nonlinearity in damp- 
ing does affect the response, especially for rounded pulse and oscillatory step inputs 
(Figures 2.7 and 2.8). Further, as more terms are included in equation (2.21), the re- 
sult progressively converges to the numerically integrated solution. It is evident from 
Figures 2.6 and 2.7, that first order solution is. good enough for case (a), but second 
order solution is needed for case (b) to match well with the numerical solution. Similar 
results were obtained over a large range of 7 up to 100 . On the other hand, for case(c) 
as seen from Figure 2.8 the two terms solution is still far away from the numerically 
integrated results. While varying 7 for this case, it was observed that only when the 
value of 7 is around 10 the closed form solution matched closely with the numerical 
solution. For a value of 7 around 10 the effect of nonlinear damping is not very much 
pronounced. The closed form solution again does not yield good result if the severity 
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Figure 2.5: Comparison of analytical and numerical results for the cubic nonlinearity 
in stiffness only. £ = 0.1, 7 = 20, e = 0.01, <7 = 1, d>i = 0. O zeroth order, * firs 

order, > second order, numerical. Curve (a): the rounded stej 

displacement; curve (b): the rounded pulse displacement; curve (c): the oscillator 
displacement step. 
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Figure 2.6: Comparison of analytical and numerical results for the cubic nonlinearity i 
damping only for the rounded step displacement. £ = 0.1, e = 0.01, =0, = 

'Y 20 O zeroth order, * first order, > second order, 


numerical 




Figure 2.7: Comparison of analytical and numerical results for the cubic nonlinearity i: 
damping only for the rounded pulse displacement. C = 0.1, e = 0.01, e\ =0, 5\ = 
7 = 20 O zeroth order, * first order, t> second order, numerical. 




Figure 2.8: Comparison of analytical and numerical results for the cubic nonlinearity 

damping only for the oscillatory displacement step. £ = 0 . 1 , € — 0.01 Cl - q S - 

7-20 O zeroth order, * first order, > second order, 


numerical 




Figure 2.9: Comparison of analytical and numerical results for the cubic nonlinearity in 
damping only for the oscillatory displacement step. C = 0- 1 5 e = 0.01, = 0, C = 1 
7 = 1 O zeroth order, * first order, t> second order, numerical. 





parameter 7 is around 1; in fact, as shown in Figure 2.9, the solution diverges as more 
number of terms are taken in the closed form solution. This is expected near 7 = 1 
with the oscillatory input (see equations (2.10) and (2.22)), because the oscillations in 
the input are at the linear natural frequency of the isolator which in turn gives rise to 
secular terms in the perturbed solutions (see equations (2.23), (2.30) and (2.34)). 

2.5.2 Effect of Nonlinear Damping 

The nonlinear damping term makes the isolator response differ appreciably from the 
response of a linearly damped isolator (e.g., see Figure 2.6). This difference increases as 
the severity parameter 7 is increased. The effect of nonlinear damping on the response 
(relative displacement) to the rounded step input is highlighted in Figure 2.10. 

These results for various values of 5* have been obtained through numerical inte- 
gration of equation (2.22) with e* = 0. With increasing values of 8 *, the peak response 
values increase up to 8* « 0.05, thereafter an increase in 5* reduces the peak value. 
However, the peak values are always higher for 8* > 0 than that for 5* = 0. Similar 
trends were also observed for other types of inputs, namely, inputs of cases (b) and 
(c). In fact this trend was also observed for any power-law damping forces with the 
exponent greater than one. 

Other Performance Variables 

The relative displacement of the mass is only one of the measures of the isolator 
performance and it is important when the rattle space is the design criteria. Other 
performance criteria of a shock isolator are indicated by the absolute displacement 
(a^/Timax), the relative velocity (A'), the absolute velocity (x' 2 /x\ max ) and acceleration 
[x' 2 /ximax) of the mass. The absolute displacement exhibits features similar to those 
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shown in Figure 2.10 ( for the relative displacement ) and is not shown again. The 
effect of nonlinear damping on the other three criteria are depicted in Figures 2.11, 
2.12 and 2.13. higure 2.11 shows that the peak value of the relative velocity is reduced 
by increasing the coefficient of the nonlinear, cubic damping term. However, Figures 
2.12 and 2.13 clearly indicate that an increase in the same coefficient results in higher 
peaks in both the absolute velocity and the acceleration of the isolated mass. 

With a shock displacement to the base, the initial velocity across the damping 
element is quite large. The cubic nonlinearity in damping then causes a large excitation 
to the mass which in turn accounts for the high value of its response. 

Quadratic Damping 

Fluid dampers, in certain range of velocity across them, can be better modelled by 
including a quadratic damping term [75]. Some elastomeric isolators are also modelled 
by a combination of linear and quadratic damping forces. In the first case the coefficient 
of the quadratic term is positive, whereas in the latter, this coefficient is a small 
(compared to that of the linear term) negative number. Both these types of damping 
can be incorporated in equation (2.2) if the term 5* (A') 3 is replaced by 5* | A' | A' 
with 5* positive or negative as the case may be. Then the equation of motion for the 
quadratic damping may be written as 

A" + 2(A' + <5* | A' | A' + A + e*A 3 = f(T) (2.73) 

Equation (2.73) was numerically integrated for various combination of parameters. It 
was seen that the effects of such quadratic damping with a negative 5* are negligible 
so far as the relative and absolute displacements are concerned. However, the absolute 
acceleration do get effected by the quadratic damping. Figure 2.14 shows that the peak 
value of the absolute acceleration can be very effectively controlled with a negative 5* . 
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Chapter 3 


RESPONSE OF NONLINEAR 
DISSIPATIVE SHOCK 
ISOLATORS TO FORCE 
EXCITATION 

3.1 Introduction 

In the previous chapter the effect of nonlinear damping on shock isolators has been 
studied for base excitation. In this chapter the effect of nonlinear cubic damping on 
shock isolators is presented for force excitation. Two types of force excitations, namely 
rounded step and rounded pulse, have been considered. Any sudden increment of the 
force may be considered as a rounded step excitation whereas a sharp rise followed by 
a gradual decline can be modelled as a rounded pulse excitation. The anvil of a drop 
forging machine is subjected to the second type of excitation. 

Unlike for the base excitation of nonlinear shock isolators, it was found that for the 
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f(t) 


x(t) 



Figure 3.1: Non-linear force excited system 

force excitation, the nonlinear damping does not deteriorate the isolator performance 
for both rounded step and rounded pulse force excitations. 

3.2 Equations of Motion 

A single degree-of-freedom shock isolator system is shown in Figure 3.1, where the mass 
is subjected to a shock force. The spring is considered as linear and the the damper 
of the isolator is taken to be nonlinear with a cubic nonlinearity superimposed on a 
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( 3 . 1 ) 


linear term. The equation of motion for the mass m, to be isolated, is 

mx + c Q x + cii 3 + k 0 x = f{t) 

where x is the absolute displacement of the mass m, and the dots denote derivatives 
with respect to time f; k 0 is the linear spring constant. c 0 and c x are the linear and 
nonlinear damping coefficients and f(t) is the shock force excitation. The initial con- 
ditions on x are taken to be zero, i.e., x — x = 0 at t = 0. Equation (3.1) is written in 
the non-dimensional form as 


X" + 2CA" + 5*(X'Y + X = F(T) 


( 3 . 2 ) 


where, X = x/x max is the non-dimensionalised displacement of the mass m; 
x max = f m ax/kch fmax ks the maximum shock force applied to the mass m, the non- 
dimensionalised time T = uiQt with w 0 = the primes denote derivatives with 

respect to T, and the other non-dimensional parameters are 


/■ Cp T* — Cj UJo ('Xmax ) * nnrl 

^ ~~ 2mw 0 } d ~~ m ’ dnC 


F{T) = 


m 

fmax 


( 3 . 3 ) 


Types of Input Shock Forces 

Two types of shock forces, rounded step and rounded pulse, have been considered. 
The functional form of these forces are similar to the functional form of the base 
excitations discussed in the previous chapter. These functional forms are continuous 
and differentiable at the time origin. The rounded step shock force represents any 
sudden increment in the force applied to the isolator mass m. The rounded pulse 
shock force represents a sudden increment in the force followed by a gradual decline. 
These forces are described mathematically as given below: 
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Case (a) Rounded step shock force: 

This type of force excitation, shown in Figure 3.2, is expressed as 

m = /m a x[l-(l + 7^)e (_7Wot) ] for t > 0 (3.4) 

with 7 as the severity parameter. 

Case (b) Unidirectional rounded pulse shock force: 

This type of force excitation is shown in Figure 3.3 and is written as 

/(t) = /™«.(e 2 /4)(7Wo() 2 e < -''"” 1) for t > 0. (3.5) 


Solution Procedure 

Equation (3.2) can be solved by using the perturbation method explained in the pre- 
vious chapter. .The non-homogeneous term F(T) in equation (3.2) may be expressed 

in generalized form using equations (3.3), (3.4) and (3.5) as given below: 

r n i 


F(T) = C + e aT 


£ A mT m 


(3.6) 


Lra=0 J 

where, C = 1, a = —7, N — 1, A 0 = — 1 and A\ = —7 for rounded step and 
C — 0, a — —7, N — 2, ^0 = 0, Ai = 0 and A 2 = e 2 7 2 /4 for rounded pulse. 

3.3 Results and Discussions 


The effect of nonlinear damping for rounded step and rounded pulse types of shock 
forces are discussed below. The linear damping coefficient £ = 0.1 and the nonlinear 
damping coefficient 5* — 0.05 have been considered. It was observed that the effect 
of nonlinear damping is discernible at high severity parameters for rounded step force 
excitation. But, for the rounded pulse force excitation, the nonlinear damping effect 
is appreciable at low severity parameters, particularly at 7 = 1. So, 7 = 50 for the 
rounded step excitation and the 7 = 1 for the rounded pulse force excitation have been 
used in all the results. 
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3.3.1 Effect of Nonlinear Damping 

Rounded Step Force Excitation 

The effect of nonlinear damping on displacement, velocity and acceleration for 
rounded step force excitation are shown in Figures 3.4, 3.5 and 3.6 respectively. These 
results are obtained by numerical integration. Unlike in rounded step displacement 
base excitation (see Figures 2.9-2.13), nonlinear damping in rounded step force exci- 
tation reduces displacement, velocity and accelerations peaks simultaneously. Similar 
to rounded step base excitation, nonlinear damping effect is appreciable when severity 
parameter is very high. At low severity parameters, the effect of nonlinear damping 
for rounded step force excitation is negligible and it is not shown here. 

Rounded Pulse Force Excitation 

The effect of nonlinear damping on displacement, velocity and acceleration for 
rounded pulse force excitation are presented in Figures 3.7, 3.8 and 3.9 respective- 
ly. Unlike in rounded pulse displacement base excitation, nonlinear damping effect on 
various response variables is insignificant at high severity parameters for rounded pulse 
force excitation. But, at low severity parameters, the nonlinear damping effect is appre- 
ciable on various response variables for rounded pulse force excitation. The nonlinear 
damping reduces the peaks of displacement, velocity and acceleration simultaneously 
for rounded pulse force excitation at low severity parameter (7 = 1). 

In overall, the presence of nonlinearity in damping does not deteriorate the per- 
formance of isolator in shock force excitation. Assumption of linearity in shock force 
excitation gives conservative estimate of the isolator performance, which is contrary 
to the base shock displacement excitation. So, in forthcoming chapters, force shock 
excitation is not considered. 

ENTRAL UB8A1V 

i. 1. ?., KANFUQ ' 

131994 
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Chapter 4 


PERFORMANCE OF 
NONLINEAR ISOLATORS AND 
ABSORBERS TO SHOCK 
EXCITATIONS 

4.1 Introduction 

Insertion of a vibration isolator between the source and receiver of vibration, and 
an attachment of a vibration absorber (a secondary system) to the main vibratory 
system are two very common methods of vibration control [1-3]. Under steady-state 
excitations, both the isolator and absorber are modelled as linear systems. However, 
for a shock loading on the system, the nonlinear characteristics present in the stiffness 
and damping properties of these systems, viz., air springs, elastomeric dampers and 
wire-rope isolators [4], come into play. In the second chapter, a single-stage shock 
isolator comprising of a parallel combination of a nonlinear spring and a nonlinear 
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damper has been considered. It was shown that the nonlinearity in the damping rather 
than that in the stiffness has more pronounced effect on different indices expressing 
the performance of a shock isolator. The presence of a nonlinear velocity dependent 
damping term with a positive coefficient was found to have detrimental effects on the 
isolator performance. 

In this chapter, different alternatives are considered for improving the performance 
of an isolator having a nonlinear cubic damping over and above the usual viscous 
damping. Since the nonlinearity in the isolator restoring force does not have any 
appreciable effect in its performance [32], the stiffness of the isolator is assumed to be 
constant in this chapter. The following four different alternatives are being considered 
and compared [77] : 

(i) inclusion of a Coulomb damper in parallel to the original nonlinear isolator, 

(ii) a three element isolator where the nonlinear damper is elastically (rather than 
rigidly) connected, 

(iii) attachment of a secondary system through a linear spring and a nonlinear 
damper, 

(iv) a two stage isolator, where each stage comprises of a linear spring and a non- 
linear damper. 

Three different types of shock inputs are considered as the base motion to be iso- 
lated. Three different indices are used to judge the overall performance characteristics 
of the isolator. Overall, it is seen that the three-element and two-stage isolators are 
preferable in the presence of nonlinear cubic damping. 
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4.2 Equations of Motion 

In this section, equations of motion for the four systems described above are derived 
in non-dimensionalized form. 


4.2.1 Type 1 Isolator : Shock Isolator with a Coulomb Damper 

Let a rigid mass m be isolated from the base excitation, x x (t), by using three elements 
in parallel viz., a linear spring, a nonlinear velocity-dependent damper and a Coulomb 
damper as shown in Figure 4.1. The equation of motion for the mass m, is 

mx 2 + c 0 (i ;2 - x x ) + c x (x 2 - iq) 3 + c f sgn(x 2 - iq) + k 0 (x 2 - xq) = 0 (4.1) 


where x x and x 2 are the absolute displacements of the base and the mass in, respectively, 
and the dots denote derivatives with respect to time t. cj is the coulomb damper 
coefficient. The initial conditions on x 2 are taken to be zero, i.e. , x 2 = x 2 = 0 at t = 0. 
Equation (4.1) can be written in the non-dimensional form 

A" + 2(A' + 5*(A') 3 + 5 f sgn( A') + A = f(T) (4.2) 


where, A = (x 2 — xi)/x x max is the non-dimensionalized relative displacement between 
the mass m and the base; the non-dimensionalized time T = uot with lu 0 = the 
primes denote derivatives with respect to T, and the other non-dimensional parameters 
i* = aatW, i, = and f(T) = -^~ 


milmatW o ' 


T 1 max 


4.2.2 Type 2 Isolator : Three Element Shock Isolator 

A three element mounting comprised of two springs and a nonlinear damper is shown 
in Figure 4.2, where an elastically connected nonlinear damper is put in parallel to 
another spring. 
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The equation of motion for the rigid mass m may be written as 


mx 2 + k 0 (x 2 - Xi) + k 2 (x 2 - x 3 ) =0. (4.3) 

Further, because the dashpot and the spring experience the same transient force, 

k 2 (x 3 - x 2 ) = co(£i - x 3 ) + Ci(x i - x 3 ) 3 . (4.4) 

The initial conditions on x 2 and x 3 are taken to be zero, i.e., x 2 = x 3 = x 2 = x 3 = 0 
at t = 0. Equations (4.3) and (4.4) are written in the non-dimensional forms as 

A" + A + n( A -Z) = f(T) (4.5) 

and 

n(Z - A) + 2(Z' + 5*Z ,z =0 (4.6) 

where, Z = (x 3 — X\)/x imax , n — k 2 /ko, and the other parameters are already defined 
in section 4.2.1. 

4.2.3 Type 3 Isolator : Transient Vibration Absorber 

Let, a secondary system be connected to the primary mass m - which is being isolated 
from the base subjected to a shock excitation. The secondary system (transient vibra- 
tion absorber) consists of a rigid mass m 2 connected to m, through a linear spring and 
a nonlinear clamper as shown in Figure 4.3. The equations of motion for the masses m 
and m 2 are 

mx 2 + co(.i '2 ~ x i ) + ci (x 2 — .x'i) ,! + ko{x 2 — aq) + C 2 (x 2 — x 3 ) 

+c 3 (i;2 — X 3) 3 + k 2 {x 2 — x 3 ) = 0 (4.7) 

and 

rn 2 x 3 + C ‘2 (:i; 3 — x 2 ) + c 3 (x 3 — X 2) 3 + k 2 (x 3 — x 2 ) = 0 (4-8) 
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where, x u x 2 and x 3 are the absolute displacements of the base, the isolated mass m 
and the absorber mass m 2 , respectively. The initial conditions on x 2 and x 3 are taken 
to be zero, i.e., x 2 = x 3 = x 2 = x 3 = 0 at t = 0. Equations (4.7) and (4.8) may be 
written in the non-dimensional forms as 


A" + 2CA' + 6*( A') 3 + A - [nY + 2( 2y /njIY' + huyfU/riY*] = f(T) (4.9) 


and 

Z" + (n/ //,) Y + 2Ca y^/tiY' + 5 2 yfvfaY ,Z = f(T) (4.10) 


where, Y = (x 3 - x 2 )/x lmax , //, = m 2 /m, u 2 = y4; 2 /m 2 , C2 = 5^, & 
and the other parameters have already been defined. 


^ 3^2 {%lmax 

m2 


4.2.4 Type 4 Isolator : Two Stage Shock Isolator 

A two stage isolator system is shown in Figure 4.4, where the main mass m (to be 
isolated) is connected to an intermediate mass m 2 (called the first stage), which, in 
turn, is then connected to the base (referred to as the second stage). In both the 
primary and secondary stages the spring elements are linear whereas the dampers are 
nonlinear. The equations of motion for the mass m and m 2 are 

'inx 2 + Co (:i'- 2 ~ -hi) 4- Ci(x 2 — i'3) 3 + ko(x 2 ~ £3) = 0 (4-11) 

and 

m 2 .'i;;j + Co (x 3 — x 2 ) + cq (2; 3 — x 2 )' 5 + ko(x 3 — x 2 ) + c 2 (x 3 — iq) 

+c 3 (x 3 — X1) 3 + & 2 (x 3 — £1) = 0 (4-12) 

where, X\, x 2 and x :i are the absolute displacements of the base, the mass m and the 
mass rn 2 , respectively. The initial conditions on x 2 and x 3 are taken to be zero, i.e., 
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X 2 = x 3 = x 2 - x a - 0 at t = 0. Using the non-dimensional parameters defined earlier, 
equations (4.11) and (4.12) may be written as 

A" - [y + 2C r + S-Y' 3 } = f(T) (4.13) 

and 

Z" + (n/ f i,)Z 4- 2C 2 \fnhiZ' + b 2 ^[^faZ' 3 + [Y + 2 (Y' + 5*Y ,3 ]/n = f(T). (4.14) 

Types of Shock Displacements of The Base 

The three commonly considered shock displacements of the base, taken as input xi (t), 
for all the above four types of isolators are already discussed in chapter 2. 


4.3 Results and Discussions 

In this section, the performance of each of the four types of shock isolator (discussed 
in previous section) systems are evaluated using the following three indices: 

Relative Displacement Ratio (RDR) = | A \ max ; 

Shock Displacement Ratio (SDR) = 

Shock Acceleration Ratio (SAR) = \ x j } max . 

The last index is expressed in dB scale by rewriting SAR = 20 log X0 ^^ max - . 

These indices are obtained for each type of isolator by numerically integrating the 
corresponding equation of motion. The Runge-Kutta Cash-Karp routine [78] has been 
used for numerical integration. In all the cases considered for the original isolator, 
the linear damping coefficient is taken as ( = 0.1 and the nonlinear cubic damping 
coefficient is taken as 5* = 0.01. 
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Figure 4.5: SDR and RDR curves for type 1 isolator as a function of 7 for the rounded 
step displacement. I. linear damping ( = 0.1, II. linear + cubic damping ( = 0.1,5* = 
0.01, III. linear + Coulomb damping ( ~ 0 . 1 , 6 / = 0.5 and IV. linear + cubic -f 
Coulomb damping < = 0.1,5* = 0.01,5; = 0.5 (suffix S for SDR and R for RDR). 
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Figure 4.G: SAR curves in dB for type 1 isolator as a function of 7 for the rounded step 
displacement. I. linear damping ( = 0.1, II. linear + cubic damping ( = 0.1,5* = 0.01, 
III. linear + Coulomb damping ( = 0.1,5/ = 0.5 and IV. linear cubic + Coulomb 
damping ( = 0.1, 5* = 0.01, Sr — 0.5 




4.3.1 Type 1 Isolator 

In this section, the effect of the Coulomb damping on the performance of a shock 
isolator system (equation (4.2)) is discussed. For the rounded step base excitation 
given by equation (2.4), the SDR, RDR and SAR curves are shown in Figures 4.5 and 
4.6, as a function of the severity parameter 7. 

It is evident by comparing curves I (linear damping) and II (linear + cubic damping) 
in Figures 4.5 and 4.6, that even for a reasonable severity parameter (7 > 5), the cubic 
damping is always detrimental to the isolator system. If a Coulomb damper is added 
in parallel to the linear viscous damper, then both the SDR and RDR arc reduced as 
is evident by comparing the cases I and III in Figure 4.5. The SAR curve in Figure 4.6 
shows somewhat worse response for low 7, i.e., the SAR. values (case III) are high for 
small 7 in the presence of Coulomb damping. This is due to the combination of a high 
value of the Coulomb damping coefficient (5/ = 0.5) and a low excitation level for small 
7. With the inclusion of a Coulomb damper to a (linear + cubic) damped isolator, the 
SDR and RDR values (curve IV) asymptotically approach those without the Coulomb 
damper from below with increasing 7 (Figure 4.5). However, in the case of SAR, the 
values without the Coulomb damper are approached from above when the Coulomb 
damper is incorporated (Figure 4.6). This is due to the high relative velocities at high 
severity parameters, which give rise to a high cubic damping force as compared to the 
constant Coulomb damping force. 

In Figure 4.5, for the cases I, II and IV, SDR « 1 + RDR, but for case III, SDR 
~ RDR with increasing 7. This is due to the fact that in the presence of Coulomb 
damping, the system comes to a rest before completing the first half cycle and the 
peaks of the relative and absolute displacements are approximately equal. 

It is obvious from the above discussion, that if SDR is the design criterion for 
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a shock isolator system, then the inclusion of a Coulomb damper is helpful. It was 
observed that il the severity parameter is indeed high, then increasing the Coulomb 
damping reduces the SDR values. For example, for 7 = 50, the parameter combination 
C = 0.1,5* = 0-01,5/ = 0.5 yields SDR = 1.994, whereas with C = 0.1,5* = 0.01,5/ = 
1.0 the SDR value is reduced to 1.598. 

The response results with the other two types of inputs, given by equations (2.6) 
and (2.7) are qualitatively similar to what has been just discussed and therefore, are 
not repeated. 

The normal stick-slip phenomenon in the presence of Coulomb damping was not 
observed with those transient excitations. Once the system sticks, it does not slip 
again. 


4.3.2 Type 2 Isolator 

Snowdon [4] has studied a three element isolator with linear elements and concluded 
that its performance is hotter than that of a two element isolator. This is specially so 
at high values of the severity parameter. In this section, the performance of a three 
element isolator in the presence of cubic damping is discussed. Figures 4.7 and 4.8 
show the SDR, the RDR and the SAR curves as a function of the severity parameter 
7 for the rounded step excitation. In these Figures, the curve V refers to the three 
element isolator, for n = 1, with all other parameter values same as in the previous 
section. It may be noted that, a three element isolator reduces the SDR and the RDR 
values close to those with the linear two-element isolator, even in the presence of cubic 
damping. This is true for all values of the severity parameter. 

Figure 4.8 reveals that for 7 > 10, the SAR values with the three element isolator 
are much lower than those with the linear two element isolator. The effect of the spring- 
ratio parameter, n, on the SDR is shown in Figure 4.9 for 7 = 50. It is clear from 
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Figure 4.7:/sDR and RDR curves for type 2 isolator/as a function of 7 for the rounded 

fa- f yft l MXtJ-rS 

step displacement. I. linear damping C = O.R II. linear + cubic damping < = 0.1,<$* = 

■K ty/lx 1 wfUJw A 

0.01^ V. type 2 isolator C = 0.1, 6* = 0.01, n = 1 (suffix S for SDR and'R for RDR). 









this Figure, that the optimum value of n is near 1. Furthermore, the value of n can 
be somewhat laigei than this ii-opUmum without much appreciable effect but reducing n 
from 7 1 optimum is undesirable. 

Comparison of Figures 4.5 and 4.6 with Figures 4.7 and 4.8, respectively, reveals 
that, at high severity parameters, the performance of a three element isolator is better 
than that of the type 1 isolator. Therefore, to nullify the cubic damping effect, it is 
better' to mount the cubic damper elastically for a single-stage shock isolator. 

Similar trends were observed for all types of excitations. 

4.3.3 Type 3 Isolator 

Snowdon [4] has also studied a, transient vibration absorber attached to an isolator 
with linear elements and concluded that when the absorber is of optimum design, the 
motion of the mounted item decays very rapidly with time. He also concluded that the 
maximum acceleration of the mounted item is considerably reduced by the dynamic 
absorber when 7 is large. In this section, the effects of attaching a vibration absorber 
over and above the nonlinear isolator are discussed. Figures 4.10 to 4.13 show the SDR 
and the SAR curves as a function of the severity parameter 7 for different types of 
excitation. In these Figures, curve VI corresponds to a transient vibration absorber, 
for //, = 0.5, C2 = 0.1 and <V = 0.01, with all other parameter values same as in the 
previous section. 

As shown in Figure 4.10, the vibration absorber docs reduce the SDR values for 
a rounded step displacement excitation (compared to that for the linear plus cubic 
damped isolator) at high values of the severity parameter (7 > 10). This reduction is 
large at high severity parameters, but the SDR values are still much higher than those 
with the linearly damped isolator. However, for the oscillatory step excitation (Figure 
4.11)), the SDR values are lower than even that with the linearly damped isolator at 
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Figure 4.10:y(SDR curves for type 3 isolator/as a function of 7 for the rounded step 

(r\. tyf<l /Ml IaJ'W 

displacement. I. linear damping ( = O.ljf II. linear + cubic damping ( = 0.1,5* = 0.0# 
VI. type 3 isolator C = 0.1, 6 * = 0.01, n = 1, S 2 = 0.01, p, = 0.5 




Figure 4.11:/(SDR curves for type 3 isolator^ a function of 7 for the oscillatory 
displacement step. 1. linear damping C = 0.}/ hnear + cubic dam P in § C “ 0-1,5 ~ 

jjtV iW 1 l ^ ' 

0.0V VI. type 3 isolator < = 0.1, <5* = 0.01, n = 1 ,62 = 0.01, n - 0.5 
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Figure 4.12:/SAR curves in dB for type 3 isolators a function of 7 for the rounded step 

M- ty*- 1 a*, i 

displacement. I. linear damping C = 0.1,/, II. linear 4- cubic damping £ = 0.1,5* = 0.01<f w'ZJw 
VI. type 3 isolator C = 0.1, 5 * = 0.01, n = 1, S 2 = 0.01, ji = 0.5 
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Figure 4.13:/SAR curves in dB for type 3 isolator's a function of 7 for the rounded pulse 

A yf^i CrJlc^-W _ ^ fattyfKl 

displacement. I. linear damping ( = 0.1/ II- linear + cubic damping C = 0.1, S = O.Oy ,^7^ 


VI. type 3 isolator C = 0.1,5* = 0.01, n = M2 = 0.01, m = 0.5 




Figure 4.14: Curve (a): The effect of stiffness ratio (n) on the SDR in type 3 isolator, 
for the rounded step displacement at 7 = 50 (( = 0.1,5* = 0.01, 5 2 = 0.01, /i = 0.5). 
Curve (b): The effect of mass ratio (jj.) on the SDR in type 3 isolator, for the rounded 
step displacement at 7 = 50 (C = 0.1,5* = 0.01, 5 2 = 0 . 01 , n = 1 . 0 ) 




all severity parameters (7 > 0.5). From the results obtained (not presented here) the 
SDR curve for the rounded pulse is similar to that with the rounded step and the RDR 
curves for all the three types of excitations are qualitatively similar to the SDR curves. 
The SAR curves for the rounded step and the rounded pulse are shown in Figures 4.12 
and 4.13, respectively. From these Figures we may conclude that the SAR values are 
not reduced by using a vibration absorber. Thus, inclusion of a vibration absorber may 
be useful only for the oscillatory step excitation, when the SDR and RDR values are 
reduced. 

The effects of the mass ratio (//,) and the stiffness ratio (n) on the SDR are shown 
in Figure 4.14. If is readily seen that for //. > 0.7, the SDR is almost independent of 
//,; whereas the SDR increases sharply with decreasing ft, for /i < 0.7. So far as the 
stiffness ratio is concerned, tilt; SDR is minimum for n = 0.8. It was observed that this 
optimum value of n remains close to 0.8 for all values of /i. The RDR and the SAR 
values were found f.o lie insensitive to the changes in the mass and stiffness ratios. 

From a. parametric investigation (not presented here), it has been found that the 
effects of the other secondary system parameters ((2 and 62) on the SDR, RDR and 
SAR are not appreciable. The values of the SDR, RDR and SAR do not depend on 
whether the secondary system is linear or nonlinear. Thus, the absorber system may 
as well be linear. 

4.3.4 Type 4 Isolator 

In this section, the performance of a two stage isolator is compared with that of the 
(linear- + cubic) damped isolator. Figures 4.15 and 4.16 show the SDR and SAR 
curves, respectively, as a function of the severity parameter 7. In these Figures, curve 
VII refers to a two stage isolator, for ft, = 0.1, with all other parameter values same as 
in the previous section. 
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As shown in Figure 4.15, t;hc SDR curve for a two stage isolator is close to that of 
the linear isolator at all severity parameters for the rounded step excitation. The SAR 
curve for the two stage isolator is lower than that of the linear isolator at all severity 
parameters for the rounded step excitation as shown in Figure 4.16. The RDR curve 
(not shown here) is similar to the SDR curve for the rounded step excitation. For the 
rounded pulse and the oscillatory step excitations, the SDR, RDR and SAR curves (not 
shown here) were observed to be similar to those with the rounded step excitation. In 
the presence of cubic damping in the primary system, a two stage isolator is better 
than all other systems for reducing the SDR, RDR and SAR. simultaneously at high 
severity parameters for all the three types of base excitations. 

The effect of the mass ratio (//.) on the SDR is shown in Figure 4.17. For m < 0.6, 
the SDR is rather insensitive to the value of but for /i > 0.6, the SDR increases 
linearly with increasing mass ratio. This effect is opposite to what has been observed 
in the type 5 isolator. 

It was further noted from the results not included here is that the SDR reduces 
with G a.n<I n. The effect of the nonlinear damping of the secondary stage on SDR 
is negligible. Thus, like in the previous case, the secondary system may be taken as 
linear. 
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Chapter 5 


2-DOF HEAVE-PITCH MODEL 
OF A VEHICLE FOR SHOCK 
EXCITATION 

5.1 Introduction 

In the previous chapters, we studied the effect of nonlinear damping on shock isolators 
in a single (heave) mode of vibration of a quarter-car vehicle model. In continuing with 
the investigation in a systematic manner, the next logical step would be to consider 
a two-dimensional (2-D) or half-car model, which includes vehicle pitch in addition 
to the sprung mass heave mode. The details of this 2-D model are given in the next 
section. 

In this chapter, two of the four isolator (suspension) models discussed in chapter 4 
are considered. The first one is a suspension system consisting of a two element isolator 
and the second one is elastically coupled three element isolator. With some special 
choice of parameters, equations are reduced to simpler forms amenable to analytical 



solutions. For general nonlinear, coupled equations, numerical results are obtained with 
the rounded displacement step is considered as the input from the road bump. Like 
in a single (heave) mode SDOF model, it was observed that in this Heave-Pitch model 
also, nonlinear symmetric dissipative damping with positive coefficient is detrimental 
to shock isolation. The elastically coupled three element isolator nullifies the adverse 
effect of the nonlinear damping in this Heave-Pitch model as well. 

5.2 System Models 

5.2.1 Model I 

The 2-DOF (Heave-Pitch) vehicle model with two element isolators is shown in Figure 

5.1. The body structure of the vehicle is represented by a rigid body of mass m with 
moment of inertia I with respect to the centroidal axis perpendicular to the vehicle’s 
plane of symmetry. The distances from the center of mass, to the front and rear axles, 
measured in the horizontal direction are a and 6, respectively. The variables y\ and y 2 
denote base inputs due to road irregularities at the front and rear tyres, _ respectively. 
It is assumed that the rear tyres travel over the same path as the front tyres; hence y 2 
is a delayed version of jq, i.e., y 2 = V\{t — to), where to = 9 y^, V is the forward speed, 
and a + b = L is the vehicle wheel base. The body motion is described by the vertical 
displacement x of the center of mass and the angle of rotation 6 about the centroidal 
axis, from the static equilibrium position shown in Figure 5.1. 

From the previous chapters, we know that the effect of nonlinearity in stiffness on 
shock isolation performance is not predominant, so, here in model I shown in Figure 

5.1, springs (/ci and k 2 ) are considered linear. The dampers, however, are considered 
nonlinear with cubic nonlinearity superimposed on a linear term. Let X\ and x?, be 
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Figure 5.1: Heave-Pitch Model I 






the vertical upward displacements, from the equilibrium position, at the points 1 and 
2 respectively in Figure 5.1. The damping forces in the front and rear suspensions are 
fdi = chjAi + cn Ai 3 and f d2 = c 20 A 2 + c 21 A 2 3 respectively; where, A, = - y 1 
and A 2 — x 2 — y 2 are the relative displacements across the front and rear dampers in 
model I ; and dot denotes differentiation with respect to time t. It is assumed that the 
pitch motion, 6, is small. 

Equations of motion for the model shown in Figure 5.1 are 


mx + /i + J 2 = 0 

(5.1) 

and 10 + fia — f 2 b = 0 

(5.2) 


where, f\ and / 2 denote the downward front and rear suspension forces at points 1 and 
2, and are given by 

/i = k\(x\ — U\) + Cio(xi ~ V\) + c'n (xj — y/i)' ! (5.3) 

and f 2 = k 2 (x 2 - y 2 ) + c 20 (f 2 - y 2 ) + c, 2l (x 2 - y 2 )'\ (5.4) 

In addition, Xi and x 2 are related to x and 6 through 

X\ = x + ad (5.5) 

and x 2 = x - bO. (5.6) 


In order to perform a parametric study, equations of motion (5.1) - (5.6) are written in 
non-dimensionalized form. The displacement variables are non-dimensionalized using 
Umax, which denotes the maximum value of the input displacement y x at the front 


wheel. The non-dimensionalized time is taken as T — ui 0 t where, ui 0 = 

Thus the non-dimensional forms of equations (5.1) and (5.2) are 

= \/(^i + k'J/m. 

X" + F l + F 2 = 0 

(5.7) 

and R 2 Z" + Fx - F 2 - a(F x + F 2 ) = 0 

(5.8) 
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where, X 


■;Z = 


— OL 


• Ql 

Umax 1 Syntax : L 1 mL * 1 ~ x mLJ^ymax 1 ~ * mu^ymax ' 

prime denotes differentiation with respect to non-dimensionalized time T. 
Equations (5.3) and (5.4) may be written in non-dimensionalized form as 


b—a . p2 41 . rp 

1 ■ ft — mL 2 ’ ri 


—p — • p 

m^oVmax 5 




F x = a^AT - Yi) + 2Ci(X( - Y[) + 6 x (X[ - T/) 3 
and F 2 = a 2 (X 2 - V 2 ) + 2G(X' - T 2 ') + S 2 (X' - K,') 3 . 
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An 
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and F 2 
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(5.9) 

(5.10) 


From the geometry of the vehicle body, the relations between A r i, A" 2 , A r , Z and or 
may be written as 


Xi 

= X + Z(l-a) 

(5.11) 

X 2 

= X-Z( 1 + a) 

(5.12) 

X 

= [(1 + a)Xx + (1 - a)X 2 )/2 

(5.13) 

z 

= (XI — X2)/2. 

(5.14) 


5.2.2 Model II 

The 2-DOF (Heave-Pitch) vehicle model with three element isolators is shown in Figure 
5.2. In this model II, both front and rear nonlinear dampers are elastically coupled. 
The front suspension nonlinear damper is coupled with a linear spring k$ and the rear 
suspension nonlinear damper is coupled with a linear spring fc 4 . Like in model I, springs 
k ] and k 2 are considered linear and all dampers (Cio, C\ X , C 2 o, C 2X ) are considered 
nonlinear with cubic nonlinearity superimposed on a linear term. Let x x ,x 2 ,xs and x 4 
be the vertical upward displacements, from equilibrium positions, at points 1,2,3 and 
4 respectively, in Figure 5.2. The damping forces in the front and rear suspensions are 
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fdi — C10A1 + CnAj and f d2 = c 20 A 2 + c 21 A 2 3 respectively; where, Aj = x :j - y x 
and A 2 = x 4 -y 2 are the relative displacements across the front and rear dampers; and 
dot denotes differentiation with respect to time t. It is assumed that the pitch motion, 
9, is small. 

Equations of motion for the model shown in Figure 5.2 are 

mx + h+h = 0 (5.15) 

and ie + f ia -f 2 b = 0 (5.15) 

where, f x and f 2 denote the downward front and rear suspension force's at. point, s 1 and 
2 in Figure 5.2, and are given by 

/1 = ki (x { -2/1) + - x :] ) (5. 1 7) 

and j 2 = k 2 (x 2 - y 2 ) + k 4 (x 2 - x A ). (5.18) 

Equations of motion for the points 3 and 4 in Figure 5.2 may Ik; written as 

^( x 3 - £1) 4- c 10 (f 3 - 2/0 4- c u (a; 3 - 2/i) :f = 0 (5.19) 

and k 4 {x A - x 2 ) 4 - c 20 (x 4 - y 2 ) 4- c 21 (:;;,, - y 2 f () (5.20) 

After non-dimensionalization, equations (5.15) - (5.20) may be written a.s 

X" 4- Fi 4- F 2 = 0 (5.21) 

R 2 Z" 4 - Fj- F 2 - cv(F' 1 4- F 2 ) = () (5 22) 

Fl = a i(*i-*i) + niC*i-A' 3 ) (5.23) 

E2 = g 2 (A 2 — Y 2 ) 4- n 2 (X 2 — X 4 ) (5.24) 

n i(X 3 - Xi) 4- 2(1 (A)' - Y[) 4- <5i(A5( - Y() :i = q (5.25) 

and n 2 (X 4 -X 2 ) + 2( 2 (X> i -Y') + 6 2 (X>-Y^ = () (5.26) 
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where *3 = ^5 X 4 = n x = n 2 = and remaining other 

non-dimensionalized variables are same as those in model I. 


5.3 Theoretical Analysis for Special Cases 

In general, dynamic equations for both model I (equations (5.7) and (5.8)) and model II 
(equations (5.21) and (5.22)) are not amenable to theoretical analysis and also it is not 
possible to solve them analytically in closed form, because these equations are coupled 
and nonlinear. However, for some special choices of parameters, these equations are 
reduced to simpler forms. Under these conditions, some general isolation characteristics 
can be extracted without numerical integration of the equations of motion. In this 
section, we discuss these special cases. For all the dynamic equations of motion, initial 
conditions X = X' = Z = Z' = 0 at T = 0 are assumed. 

5.3.1 Symmetric System 

The simplest system is the symmetric system, in which the front and rear suspension 
characteristics are identical and the offset parameter a is zero. For R = 1 and Y 2 = Yi, 
the symmetric system behaves like an SDOF system with X = X\ = X 2 and Z = 0. 
This SDOF system is similar to the two element isolator system without nonlinear 
stiffness in chapter 2 for model I and is similar to the three element isolator system in 
chapter 4 for model II. 

From Chapter 2, we already know that the presence of nonlinear damping is detri- 
mental to SDOF two element isolator. From Chapter 4, we know that the three element 
isolator nullifies the adverse effect of nonlinear damping. 
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5.3.2 Decoupled System 

From equations (5.13), (5.14) and for 

R 2 = 1 - a 2 , (5-27) 

equations (5.21) and (5.22) may be simplified as 

(1 + <y)X" + 2F\ = 0 (5.28) 

and (1 - a) A" + 2F 2 = 0. (5.29) 

Equation (5.27) implying / = mab is called as the decoupling condition. If this 

condition is satisfied, then the 2-DOF system reduces to two SDOF systems (equations 

(5.28) and (5.29)) with mass m\ — nrJb/L at point 1 and mass m 2 = ma/L at point 2, 
which are connected by a massless rod. The above nonlinear decoupled equations (5.28) 
and (5.29) for model I can be solved analytically by using the perturbation method, 
explained in the second chapter for shock excitations. 

Without solving the decoupled equations (5.28) and (5.29); we can deduce the 
relation between X and Z for three special inputs of Y x and Y 2 . 

Case I: Y 2 = 0 for T > 0 

This case occurs in practice, when the vehicle moves forward very slowly. Equation 

(5.29) gives X 2 = 0 for Y 2 = 0 (see equation (5.9) for model I and equations (5.24) 
and (5.26) for model II). Equation (5.12) gives A" = Z( 1 + «) and equation (5.14) 
gives Z = Xi/2 for X 2 = 0. In this case it is sufficient to solve only one equation 
viz., equation (5.28), to determine the responses X and A. Heave and pitch modes are 
linearly dependent and their ratio X/Z = 1 + a is a constant. For a = 0, A r = Z\ this 
is an example of a repeated eigenvalue and eigenvectors problem in 2-DOF system and 
both natural frequencies and normal modes are identical. 
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Case II: Y 2 = Y X (T - T 0 ) for T > 0 and T 0 = Lco 0 /V > 0 

In this case, the system behaviour is similar to the previous case, up to T < To; 
it means that the input at the rear wheel Y 2 (T) = 0 for T < T 0 ; hence, X 1 2 (T) = 
0 ]X(T) = (1 + a)Z{T) and Z(T) = X x {T)/2 for T < T 0 . At T = T 0 , there will be a 
sudden change in the heave and pitch responses, due to input y 2 at the rear wheel. In 
shock excitations the maximum peak occurs during the first cycle. If T 0 is less than the 
time required to attain the maximum peak due to input y x at the front suspension, only 
then input y 2 can alter the maximum peak. So, it is important to study the response 


at low values of T 0 , which correspond to high vehicle speeds. 

Let us assume that 

Fi = (1 + a)F(Xi,Xi , Yi, Yi') (5.30) 

and F 2 = (l-a)F(X 2 ,X 2 ',Y 2 ,Y 2 '), (5.31) 

then equations (5.28) and (5.29) may be rewritten as 

X” + 2F(X l ,X l ',Y 1 ,Y 1 ') = 0 (5.32) 

and X'J, + 2F(X 2l X 2 , Y 2 , Y 2 ) = 0. . (5.33) 


Equations (5.32) and (5.33) are similar and independent of a. Since, Y 2 is a delayed 
version of Y x \ i.e., Y 2 (T) — Y X (T — T 0 ); X 2 will also be a delayed version of X x ] i.e., 
X 2 (T) = X x (T — To). 

The assumption in equations (5.30) and (5.31) is indirectly putting a condition on 
the characteristics of the front and rear suspension springs and dampers. It implies 
that, the stiffness ratios {a x /a 2 ; n x /n 2 ), the linear damping coefficients ratio (( 1 /( 2 ) and 
the nonlinear damping coefficients ratio (5 x /8 2 ) of the front to the rear suspensions are 
simultaneously equal to = 6/a. In other words, 

a x / <*2 = C 1 /C 2 = 81/82 = n x /n 2 = 6/a. (5.34) 
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This condition i.e., equation (5.34) may be called as the ‘invariant ratio’, which, 
when satisfied, renders X x and X 2 invariant to a. From equation (5.14) it can be 
deduced that, Z{T ) is also invariant to a. 

Case III: Y 2 = Y x 

This case occurs in practice when the input excitation Y\ is periodic with period 
T 0 > 0, that means both wheels are subjected to the same input simultaneously. For 
Y l - Y 2 , equations (5.32) and (5.33) are identical to each other, which implies A'j = X 2 . 
From equations (5.13) and (5.14) it can be concluded that X — X\ = X 2 and Z ~= 0. 
Even though the system is asymmetric ( b > a); in this special case Y 2 — V'j and 
under the assumption of invariant ratio given by equation (5.34), the decoupled system 
behaves as if it is symmetric ( Section 5.3.1). We know from calculus of variations that 
a symmetric system always gives optimum response. In our vehicle model, to satisfy 
the assumption of invariant ratio, front to rear stiffness and damping ratios should be 
b/a , which is always greater than unity, for forward center of mass (b > a). 

The discussions presented above for the three special cases are true for any type of 
excitation (i.e., harmonic, random and transient). 

5.3.3 Linear System 

In this section, the solution of a linear system for model I is obtained by using Laplace 
transform. This solution has been used for studying the effect of linear damping on 


vehicle suspension performance for model I. For the linear syste in (5j = S 2 = 0), 
equations of motion (5.7) - (5.10) reduce to 

X" + F\ + F 2 = 0, (5.35) 

R 2 Z" + F 1 -F 2 - a(F 1 + F 2 ) = 0, (5.30) • 

F, = a 1 (X 1 -Y 1 ) + 2( 1 (X[-Yl) (5.37) 
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and F 2 = a 2 (X 2 - Y 2 ) + 2C, 2 {X' 2 - Y'). 


(5.38) 


The solution of linear equations (5.35) and (5.36) for initial conditions X = X 1 = 
Z — Z' = 0 at T = 0 is expressed as the algebraic sum of two responses one due to 
only input \\ and the other due to only input Y 2 = Yi(T - T 0 ). 

Let Xf and Xi , be the heave mode responses due to inputs Y\ and Y 2 , respectively, 
and Zf and Z b be the pitch mode responses due to inputs Y\ and Y 2 , respectively. Now, 


the heave (X) and pitch (Z) mode responses may be expressed as 

X(T) = X f (T) + X b (T-T 0 )u(T-T 0 ) ( 5 . 39 ) 

and Z(T) = Z f (T)-Z b (T-T 0 )u(T-T 0 ), (5.40) 

where A/(T) = Xo(aiXii^nJh2, (5-41) 

A"b(T) = Xo(a 2 ,(2,huJh5,hi6,T)] (5.42) 

Z/(T) = Xo(ai,(i,h 2 i,h 22 ,h 2 3 : T)-, (5.43) 

Zb(T) = Ao(a2, ( 2 , h 24 , h 2 5, h 25 , T) (5.44) 


and u is a unit step function, and the functional form of Xq and expressions for other 
variables are given in Appendix B for input Y\ given in the next section. 

5.4 Input from The Road Irregularities 

In this chapter, a rounded displacement step has been assumed as input to the vehicle 
front and rear suspensions from the road irregularities. This rounded displacement step 
represents a sudden bump on the road. The expressions for the rounded displacement 
step in spatial co-ordinate s are given below: 

2/i (s) = Vi max [! — (!+ 7 S s/L) e (_7 ^ /L) ] for s > 0 

and y 2 (s) = y x (s - L)u{s - L) 
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(5.45) 

(5.46) 



where, 7 S is the spatial severity parameter and L is the wheel base of the vehicle. 
Equations (5.45) and (5.46) may be written in temporal co-ordinate t as, 


»,(«) = yi™«[l -(1 + 7,^/0 v ' ,/1) 

and y 2 (t) = yi{t - t 0 )u(t - t 0 ) 


for t > 0 


(5.47) 

(5.48) 


where, t = s/V and t 0 = L/V is the time delay. 
In non-dimensionalized form 


Y X (T) = [l - (1 + jsT/To) c { ~^ T/To) 
and Y 2 (T) = Yi{T - T 0 )u(T - T 0 ) 


for T > 0 


(5.49) 

(5.50) 


where T 0 = Luiq/V is the non-dimensionalized time delay. 


5.5 Results and Discussions 

The coupled nonlinear equations (5.7) and (5.8) are numerically integrated by using 
Runge-Kutta Method. The performance of the vehicle suspension is evaluated by using 
the following four indices: 

Heave Displacement Ratio (HDR) = ; 

Pitch Displacement Ratio (PDR) = ; 

Heave Acceleration Ratio (HAR) = 20 log 10 K, ; 

2^1 mar. 

Pitch Acceleration Ratio (PAR) = 20 logio >^ 5 - £ . 

^ *1 max 

In general, four wheeled vehicles have forward center of mass (a > 0) and II > 1. 
Here, we have considered a = 0.2 and R — 1.5 for the purpose of numerical integration. 
Spatial severity parameter, = 10 has been considered. The effects of linear and 
nonlinear damping on HDR, PDR, HAR and PAR are discussed. 
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5.5.1 Effect of Linear Damping 

Model I 

Numerical results are obtained for various values of damping by using the closed form 
solution given by equations (5.39) and (5.40). From the decoupled system (see section 
5.3.2), i.e., with R 2 = 1 — a 2 , if the stiffness ratio and damping ratio of front to rear 
suspensions are simultaneously equal to the invariant ratio (see equation (5.34)), then 
we get optimum X and Z. For forward center of mass (a = 0.2), the invariant ratio 
(b/a = 1^) is greater than unity; it means that high damping and high stiffness at 
the front suspension rather than at the rear suspension give the optimum response for 
the decoupled case. This is especially so at very high vehicle speed (i.e. To is close 
to zero). Figure 5.3 shows the effect of linear damping on HDR and PDR for the 
coupled system (model I). By increasing the total damping (Ci + £ 2 ) in the system, 
HDR decreases at high speeds (To < 2). For a given total damping (Ci + £ 2 ), at high 
speeds (T 0 < 3), (Ci > £ 2 ) as compared to (Ci < C 2 ) gives low HDR. At low speeds 
(T 0 > 3.5), (Ci = 0.1, C 2 = 0.2) gives lower HDR than for other damping coefficients 
(Ci = 0.1, C 2 = 0.1) and (Ci = 0.2, C 2 = 0.1). Similar observations were found for 
various values of damping coefficients. 

At high speeds (see Figure 5.3), PDR does not change much by varying damping 
either in the front or in the rear suspension. For all damping values, PDR initially 
increases with T 0 up to (To ~ 4), after that it is almost constant. At low speeds 
(To > 4.5), if the total damping (Ci + C 2 ) ia the system increases, then PDR decreases 
marginally. For a given total damping (Ci + C 2 )) reduction in PDR is more for (Ci > C 2 ) 
than with (Ci < C 2 )- 

Figure 5.4 shows the effect of linear damping on HAR for model I. At high speeds, 
IIAR increases with increasing total damping (Ci +C 2 ) i n the system. This observation 
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is similar to that for the SDOF system [4], in which SAR increases at high severity 
parameters with increasing damping value. For a given total damping (Ci +C 2 ), at very 
high speeds (To < 1), HAR is insensitive to whether (£i < £ 2 ) or (£i > £ 2 ). 

Figure 5.5 shows the effect of linear damping on PAR for model I. For a given total 
damping 4- £ 2 ), (Ci > C 2 ) gives lower PAR than that with (£i < £ 2 ) almost at all 
speeds. By increasing the total damping + £ 2 ) in the system, PAR increases with 
increasing vehicle speeds up to To < 3. 

In conclusion, for a given total damping ((\ + ( 2 ), at high speeds (Ci > C 2 ) gives 
lower HDR and PAR as compared to those with (£1 < £ 2 ). At high speeds, increasing 
the total damping in the system reduces HDR, but increases HAR and PAR. For better 
isolation, all four indices, viz., HDR, HAR, PDR and PAR should be simultaneously 
as low as possible. These are conflicting requirements which can be fulfilled only 
partially by a two-element isolator (model I). The three-element isolator model with 
linear damping is discussed next to see if both HDR and HAR can be simultaneously 
brought down. 

Model II 

In the fourth chapter (see Figure 4.7), it has been observed that the optimum SDR 
occurs at a spring ratio parameter, n = 1.0. So, here we have considered n\ = n 2 = 1.0 
in model II. Ci = C2 = 0.1 has been taken in both model I and model II. Figure 5.6 
shows comparison of models I and II for HDR and PDR. At all speeds, there is hardly 
any difference in the performance of these two models so far as HDR and PDR values 
are concerned. 

Figure 5.7 shows comparison of models I and II for HAR and PAR. At high speeds 
(To < 0.2), there is considerable reduction in both HAR (about 6 dB) and PAR (about 
10 dB) values for model II as compared those of model I. However, at lower speed 
regimes, both models yield similar values for HAR and PAR. 
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In conclusion, at high vehicle speeds model II performs better than model I by 
reducing HAR and PAR, while maintaining similar values of IIDR and FDR. 

5.5.2 Effect of Nonlinear Damping 

Model I 

From chapter two, we know that nonlinear damping is detrimental for shock isola- 
tion in SDOF system especially at high severity parameters (see section 2.4.2). Figure 
5.8 shows the effect of nonlinear damping on HDR and PDR for model I, with linear 
damping coefficients Ci = C 2 = 0-1 and nonlinear damping coefficients = 0.01. 

As in the SDOF system (see Figure 2.10), here also the nonlinear damping is detri- 
mental to the. shock isolation for model I at high speeds. It is evident from Figure 5.8 
that both HDR and PDR are increased due to nonlinear damping. At low velocities 
(To > 1), however, the effect of nonlinear damping is negligible. Figure 5.9 shows the 
effect of nonlinear damping on HAR and PAR. Both HAR and PAR increase due to 
nonlinear damping up to T 0 « 2. After T 0 > 2, again the effect of nonlinear damping 
is negligible. 

Therefore, due to nonlinear damping, all the four performance indices, HDR, PDR, 
HAR and PAR are larger at high speeds. It was observed that the distribution of damp- 
ing between the front and rear suspension does not improved the vehicle suspension 
performance unlike in linear damping case, since the presence of nonlinear damping 
itself is deteriorated the vehicle suspension performance. 

In the next section, performance of models I and II are compared in the presence 
of nonlinear damping. 

Model II 

Figure 5.10 shows the comparison of models I and II with reference to HDR and 
PDR values. For model I, Ci = C2 = 0.1 and 5 X = 5 2 = 0.01 and for model II, 
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Cl = C 2 = o.l; 5y = 5 2 = 0.01 and n x = n 2 = 1 have been used. At, high speeds (7’ 0 < 1), 
both HDR and PDR values are effectively reduced in model II when compared to those 
in model I. However, at lower speed regimes, both models 1 and II yield almost similar 
values for HDR and PDR. 

HAR and PAR values are presented in Figure 5.11 for both models I and II. At 
high vehicle speeds (To < 1.5), reduction in both HAR and PAR values is very high 
for model II when compared to those of model I. For example, this reduction is 25 
dB and 35 dB in HAR and PAR, respectively, at (T 0 = 0.2). At lower speed regimes 
(To > 1.5), the difference between models I and II for both HAR and PAR values are 
negligible. The adverse effect of nonlinear damping is completely nullified in model II. 
Like in SDOF system, it is better to mount the nonlinear dampen' elastically, in this 
2-DOF model also, to get better isolation performance. 

It is evident, by comparing Figure 5.6 with Figure 5.10 and Figure 5.7 with Fig- 
ure 5.11 respectively, that the performance of model II is similar for both linear and 
nonlinear damping cases. In conclusion, model II reduces all four performance indices 
HDR, PDR, HAR and PAR in the presence of both linear and nonlinear damping. 
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Chapter 6 


2-DOF HEAVE-PITCH MODEL 
OF A VEHICLE FOR RANDOM 
EXCITATION 

6.1 Introduction 

In the previous chapter, the performance of nonlinear suspension system has been 
studied in the context of a 2-DOF heave-pitch model of a vehicle subjected to shock 
excitation has been studied. In this chapter, the same model of a vehicle subjected 
to random excitation is investigated. The road surface profile on which the vehicle 
is traversing is modelled as the output of a first order shaping filter to white noise 
excitation. To determine the response of the vehicle, Monte Carlo simulation has 
been employed. The road surface excitation has been simulated as a seiies of cosine 
functions with weighted amplitudes at evenly spaced frequencies and random phase 
angles. Only a two element isolator model is considered. Root Mean Square (RMS) 
values of response variables are treated as performance indices. It is observed that the 
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nonlinear cubic damping reduces the RMS heave and pitch displacements but increases 
both the RMS heave and pitch accelerations. 

6.2 Equations of Motion 

The 2-DOF (Heave-Pitch) vehicle model with two element isolators, shown in Figure 
5.1, is considered again. The description of the vehicle model remains same as given 
in the previous chapter, but the input excitations y\ and y 2 are now considered as of 
stochastic nature. Equations of motion are also identical to the equations of motion 
given by equations (5.1) - (5.2). The non-dimensional forms of equations of motion 
given by equations (5.7) - (5.12) are repeated below for the sake of completeness: 


X" + F\ + F-2 = 0 


(6.1) 

R 2 Z" + Fx - F 2 - a(Fi +F 2 ) = 0 


(6.2) 

Fi = a l (X l -Y l ) + 2<; 1 (X[-Y{)+6 l (X , 1 

- n ') 3 

(6.3) 

F 2 = a 2 (X 2 - Y 2 ) + 2C 2 (X' - y 2 ') + 6 2 (X^ 

-nr . 

(6.4) 

AT = X + Z(l-a) 


( 6 . 5 ) 

AT = X - Z{ l + a). 


(6.6) 


As in the previous chapter, Y 2 {T)(= '1)2/ Umax) is considered as the delayed version 
of Yi(T){= y\/yr n ax) with Y 2 {T) = Yi(T — T 0 ), but y rnax here denotes the RMS value 
of the input displacement y v at the front wheel. All the non-dimensional variables are 
identical to those given in the previous chapter. In the next section, the statistical 
properties of y x are described. 
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6.3 Road Profile Model 


The classification of rough surfaces has received much attention over the years. 
Different models of the road surface roughness are described in the literature [79-84], 
Among various methods of analysis, the most common one is via, wave number spectra 
or covariance properties. The techniques rely on one crucial assumption that a rough 
ground profile may be regarded as the output of a white-noise excited, spatial shaping 
filter. 

Consider a road surface which is homogeneous and isotropic, with a displacement, 
profile y{s), where s is the horizontal distance, assumed to be a. member function of 
an ergodic random process {'</(«) } [80]. The profile autocorrelation function may be 
defined in terms of the horizontal lag, //,, by 

Ry{h) = E[y{s)y{s + h)) (b.7) 

where E indicates the expectation operator taken over all members of the process. 
When a vehicle travels across a road surface, the displacement experienced at its front 
wheel must be expressed as a function of time, yi(t), with magnitude equivalent to the 
profile displacement at the point of contact at that instant: thus if at time /.j the point 
of contact is at S\ we may write 

V\{h) = y(si). (G.8) 

The distance s will be a known function of t and s — Vi for a constant vehicle speed 
V. The autocorrelation function may be expressed as a function of expectations at, the 
two instants of time concerned. Thus 

R vi (*i. fe) = E[yi(h)y 2 {t 2 )] = £fo(*i)y(s 2 )] = R y {. s 1( ,s 2 ). (0.9) 

For the constant speed V, the input random process y is stationary. For a stationary 
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random process, y(s), equation (6.9) may be written as 


R yi (r)=R y (h) 


( 6 . 10 ) 


where r = t 2 - h ; h = s 2 - Si and r = h/V. 

The temporal power spectral density S yi (u>) is given by[85], 


S„M = ^ Me-“ 'dr 


( 6 . 11 ) 


and the corresponding spatial spectral density is 


Sy(») 


i r°° 


/ DO 

"OO 


cr il,h dh 


( 0 . 12 ) 


2tT j —oo 

where u = 27r/T p is the circular frequency; u — 2n/\ is the wave number; 7 ), is the 
time period and A is the wavelength. 

From equations (6.9)-(6.12), the relation between S Vl (u) and S y ( 7 ) may be written 
as 

^7/1 ( w ) = ^^(^ = u/V) (6. 13) 

The power spectral density (psd) function of road irregularity is assumed to be of the 
form [60] 


S yi H 


a 2 qV 


(6.14) 


7 r[o ; 2 + (qV) 2 ] 

where V is the vehicle forward speed, ui is circular frequency, a 2 denotes the variance 
of the road irregularity and q is a coefficient depending on the type of road surface. 
The process y\(t) described by equation (6.14) is the output of a linear first order filter 
with white noise input given by 


ih(t) + qVy i(t)=W(t). 

W(t) is a white noise process with covariance function 


(6.15) 


E{W(t)W T (t)} = 2 a 2 qV<) , - t A ) 


(6.16) 
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i x ronresf'uts tlu’ Dii'aodc'll.a iunuUon. 
pr l azotes the expectation operator and () h i - • 

W den0tGS 1 , . . n , v dvon hy equation (G.U) may he expressed m 

The psd function of road irregular ity b • Y 


terms 


of non-dimensional parameters as 


(a 2 M£ij } 

5 *( n ) = (QTo) 2 + IW 


(G.17) 


and ft ( lL. 


where, Q, - o>/u>o, # l + ** )/W equation (G.17) may be 

Similarly, the autocorrelation funct.on to, the 1-' W'"' - 


written as 


Rv 


» Vi (r) = o l exp (l -/Ir/'/o ()■ 


(G.18) 


6.4 Monte Carlo Simulation 


■ 1 r s-unnlim- experinumt.|H(),S7j. It involves a 

The Monte Carlo simulation is a statistical •« 1 t 

* f , n /| m mibsouuont (M.wminaUon ol Uu 
repeated generation of random variates a.i.1 a -l" „ 

Proper statistical processing provides esti, nates of -M »*— ’ ***^ 

The random input process V, (T) from the road «urf»* * s ‘" mW '"' * " 

nozuka’s spectral representation method[G5). 1.. this method, sample 
stochastic process are generated using a cosine series lonnuhr. 1 hose samp e 
accurately reflect the prescribed probabilistic characteristics ol the 
when the number N of the terms in the cosine series is largo. I he eusom . c av 

power spectra, density or autocorrelation function approaches hug targ . . 

(actual) function as the sample size increases. In addition, the generated sample unc- 
tions possess ergodic characteristic in the sense that the lo„.p..raU.y-avcni R c< ..ica. 
value and the autocorrelation function arc identical with the coni spending g . 
when the averaging takes place over the fundamental period of the cosine series. 
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autocorrelation function and consequently the power spectral density respectively con- 
verge at the rate 1/N 2 to the target autocorrelation and target power spectral density 
of the random process when N — > oo. The most important property of the simulated 
stochastic process is that it is asymptotically Gaussian as N — > oo by virtue of the 
central limit theorem. 

Simulation Formula 

The random input process Yy(T) with mean value equal to zero and autocorrelation 
function Ry^r) given by equation (6.18) can be simulated by the following series with 
N -> oo[65]; 

V.CT) = (v^/w) £ MV + i)| (0.19) 

71--0 

where, A n — \j2S Yi (il n )Ai). 

5 n (a0 is the psd of road surface given by equation (6.17), which is of insignificant 
magnitude (or may be assumed to be zero) below a lower limit and beyond an upper 
limit Q u . Thus, Q n is considered in the range 

4>n is the random phase angle uniformly distributed between 0 and 27 r; 

■ tt n = Qi + (n- 1/2) AO, 
iY n = il n + 5Q n , n = 1, 2, ..., N\ 

An = (Q„ - Qi)/N, 

with N being the number of intervals. A small random frequency SQ n is introduced 
to avoid the periodicity of the simulated process and is uniformly distributed between 
-Ail'/2 and AiY/2 with Ail' « AQ. The input random process Y](T) is simulated 
with N = 1000, n L = 0, n u = 50 Hz and Ail' = Ail/20. 
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The autocorrelation function oi the simulated (generated) random process )](/') is 
compared in tho next section with the actual (target) autumn elation function given 
by equation (0.18). 

6.5 Results and Discussions 

Road roughness parameter q - 0.15 and o’ 0.0 >. 10 ''nr foi asphalt road have been 
reported in several refenmees [70 8-l|. I ‘or a wheel base /, '.’.a m, p <//< value comes 

out to be 0.375. For w 0 =•• 10 rad/sec, 7’,, 1 represents vehicle speed \' 25 m/sec 

(i.e., 90 Km/h). In Figure 0.1, tin* generated autocorrelat mu function is compared 
with the target autocorrelation function given by equation (ti IK) with ,»* 0,375 and 

To = 1. It is seen that the generated autocorrelation function is matching, very well 
with the target autocorrelation function for all ;. Similar matching, was observed for 
other values of To also. Both autocorrelation functions are asymptotically approaching 
the mean value which is zero according to our assumption. 

In all the results, like in the previous chapter, t{ 1.5, n 0.9, r’ t i’ ; . 0.1 
and = 82 = 0.01 have been considered, liquations (0.1) ami (0.3) are integrated 
numerically by using Runge-Kutta method for the input, sample functions generated by 
equation (6.19). Figure 6.2 shows the effect of nonlinear cubic damping, on HMS heave 
displacement. It may be observed that nonlinear cubic dumping reduces the HMS heave 
displacement at all speeds To > 0.2. Here, high 7|) corresponds to low speeds and low / b 
corresponds to high speeds. T 0 - 0.2 represents a. speed of 125 m/see (*150 Km/h) and 
T 0 = 10 represents a speed of 2.5 m/s or, (!) Km/h) for w„ 10 rail /sec. At low speed 
regimes (To > 3), there is not much difierence between linear and nonlinear responses, 
although the nonlinear response is somewhat lower than the linear response. Both, 
linear and nonlinear responses are increased with an increase in speed up to 7jj 1.0. 
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While the nonlinear response increases monotonically with increasing speed, the linear 
response exhibits a resonance type phenomenon at 7 q 1-0. I'or I b < i, linear 
response decreases with increasing speed. This result seems to be paradoxical, since 
it is expected that the response increases with increasing vehicle speed. However, it is 
observed that beyond a critical speed, the performance of vehicle suspension improves 
rapidly with speed increment. Similar behaviour has been reported in physics literature 
[88-91] for linear and nonlinear stochastic systems, with nonlinear models consisting of 
only nonlinearity in stiffness and not in damping. This phenomenon termed ‘stochastic 
resonance’ was explained analytically only for soft Dulling oscillators having double- 
well potentials. Hitherto, there is no explanation available for this phenomenon for 
hard Duffing oscillators with a single well potential. The system under consideration 
here also has a single well potential (due to the linear springs). However, it may be 
pointed out that the critical speed for the present system, beyond which the linear 
response decreases with increasing speed, is found to be very high (beyond tin; range 
of practical interest). 

Figure 6.3 shows the effect of nonlinear damping on RMS heave acceleration. Here, 
the nonlinear response is higher than the linear response at all speed regimes. At high 
speed regimes, the difference is very high. At 7b = 1, the difference is 8 dB, but at 
To = 10 the difference is only 1 dB. 

Figure 6.4 and 6.5 show the effect of nonlinear damping on RMS pitch displacement 
and acceleration respectively, fn Figure 6.4, both linear and nonlinear responses show a 
critical speed around To = 4, which gives a peak RMS pitch displacement. The nonlin- 
ear RMS pitch displacement is always lower than the linear RMS pitch displacement. 
The difference between the linear and nonlinear responses is maximum at the critical 
speed. 

The RMS pitch acceleration (see Figure 6.5) is always more in the presence of 
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nonlinear damping. The difference between the linear and nonlinear responses is very 
high at high speed regimes. The linear response shows a critical speed at which the 
RMS acceleration reaches a maximum value, whereas the nonlinear response shows a 
monotonic increase with increasing speed. 

In summary, stochastic resonance was observed for both heave and pitch displace- 
ments wdien the damping is linear. Unlike for shock excitations, in this case the nonlin- 
ear cubic damping reduces both heave and pitch displacements. However, the presence 
of nonlinearity in the damping force increases both heave and pitch accelerations at 
high speed regimes. 
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Chapter 7 


CONCLUSIONS 

7.1 Conclusions 

In this thesis, the effect of nonlinearity in the isolator elements has been investigated 
for shock and random excitations. Both base and force excitations have been stud- 
ied. A simple analytical method has been developed using perturbation method with 
Laplace transform to determine the transient response of a single degree-of-freedom 
(SDOF) system in the presence of nonlinear, dissipative shock isolators. This method 
is applicable for any order of nonlinearity, both in the restoring and damping forces, 
expressed in the form of polynomials. However, the closed form results are obtained 
with only cubic nonlinearity both in the restoring and damping forces. Three types of 
base excitations, namely, the rounded step, the rounded pulse and the oscillatory step 
are considered. Analytical results obtained by the present method are compared with 
those obtained by direct numerical integration. Numerical results an; included for a 
typical elastomeric damper which can be modelled by a combination of linear viscous 
damping and a dissipative quadratic damping with a negative coefficient. In the pres- 
ence of inherent cubic damping in a nonlinear shock isolator with base excitation, four 
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different modifications of the isolator are considered to improve its performance. A two 
degree-of-freedom (2-DOF) model for a vehicle undergoing heave and pitch motions is 
investigated for both shock and random excitations. It is assumed that the vehicle is 
moving with a constant speed so that the rear wheels are subjected to the same input 
as the front wheels but with a time delay. Shinozuka’s spectral method is used for 
generation of the random input process. Monte-Carlo simulation has been employed 
to find the response of the vehicle. 

The major conclusions of the present thesis are listed below: 

(i) A simple analytical method has been developed using perturbation technique in 

conjunction with Laplace transform to determine the transient response of an SDOF 

¥ 

system in the presence of nonlinear, dissipative shock isolators. It is observed that 
in most cases, the closed form perturbed solution match well with the numerically 
integrated results. 

(ii) It has been established that the effect of nonlinear damping is more predominant 
than the effect of nonlinear stiffness so far as the performance of shock isolators for 
base excitation is concerned. The presence of nonlinear n th power-law damping, with 
a positive coefficient and the index n greater than unity, is harmful for base excited 
shock isolators. However, quadratic damping with a negative coefficient improves the 
performance of a shock isolator. 

(iii) With force excitation, the presence of nonlinearity in damping does not deterio- 
rate the performance of a shock isolator, rather it marginally improves the performance 
of a linear shock isolator. 

(iv) To improve the performance of a shock isolator with base excitation, in the 
presence of nonlinear cubic damping, four methods, namely, Coulomb damped isolator, 
three element isolator, vibration absorber and two stage isolator have been considered. 
Overall, a three element isolator and a two stage isolator are shown to nullify the 
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adverse effect of cubic damping nonlinearity. 

(v) 2-DOF coupled nonlinear asymmetric systems behave like a symmetric system 
for a particular choice of parameter values. At very high speeds and lor lorward center 
of mass, high damping and high stiffness at the front suspension rather than at the 
rear suspension give the optimum response for the decoupled system. 

(vi) The detrimental effect of nonlinear cubic damping in a two element isolator 
and its mitigation by using a three element isolator are exhibited by both SDOF and 
2-DOF models. 

(vii) With a random base excitation, ‘Stochastic Resonance’ type phenomenon is 
observed for both heave and pitch displacements when the damping is linear. Unlike 
for shock excitations, in this case the nonlinear cubic damping reduces heave and pitch 
displacements. However, the presence of nonlinearity in the damping force increases 
both heave and pitch accelerations. 

7.2 Scope of Future Work 

The following are the directions along which future research can be carried out: 

(i) In 2-DOF vehicle model, speed of the vehicle is assumed constant. The effects 
of variable speed with acceleration/deceleration can be looked into. 

(ii) 2-DOF heave-pitch vehicle model can be extended to heave-pitch-roll vehicle 
model with four inputs at the four wheels. Extra degrees-of- freedom can be; introduced 
to treat the sprung and unsprung masses. 

(iii) It was observed that the nonlinear asymmetric system behaves like a symmetric 
nonlinear system for a particular choice of parameter values in 2-DOF vehicle model. 
One can extend this work to investigate whether the behaviour of nonlinear asymmetric 
system (discrete or continuous) in general can be studied as a perturbation on nonlinear 
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symmetric systems. 

(iv) The phenomenon of stochastic resonance can be investigated in nonlinear con- 
tinuous systems also. 
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Appendix A 


Solution of Ordinary Second Order 
Differential Equations for Three 
Special Non-homogeneous 
Functions 


In this appendix the solution of linear second order differential equations with time- 
invariant coefficients for three types of generalized non-homogeneous functions, with 
initial conditions x = x = 0atf = 0, are presented. 

Consider the equation 

i + 2C x + x = f{t) (A-l) 

with initial conditions x(0) = x(0) = 0. 

Solution of equation (A.l) is obtained by using Laplace transformations for three 
typos of non-homogeneous functions /(£). 

Cascfi) /(t) = e“[E" = „X„t m ] 
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This is the generalized form of both rounded step and rounded pulse inputs. 
Solution of equation (A.l) may be written as 


®(f) 

where, (t) 
and X 2 (t) 

P 

a. 

i 

Qi + iP\ 

B N+i 

Bn 

Bn+ 1— m 


= e 


*1 (i) +x 2 (t) 

^ [i°i cos(]0t) + Qi sin(/3i)] / (3 


= e 


,a£ 


E B m+1 r/m! 

_77l=0 


-C + *0; 


-i); 


/VI 1 . 


y ]{N - m)\AN~,ni a a )" 1 / (°' “ a ) 

.m=0 

= N'.Ajv/T; 

= [(AT - l)!A Ar _i -2(a + C)^N+i]/f\ 

= [(A — m)!Aiv-m ~ 2(a 4- m ~ Bn+x- m] /' 1 


a 2 + 


for 2 < 7//. < N. 


Case (ii) f{t) = e at {C m cos(w m t) + D m sin(cj m i)}J 

This function is a generalization of the oscillatory step input. Solution ol equation 
(A.l) then may be written as 


x{t) 

where, X\ (t) 
and x 2 (t) 


= e 


X\{t) + x 2 {t) 
e ~ t ’ t [P 2 cos (/it) + Q‘i sin 

N 

y cos(w m t) + 5 m sin(w m t) } 

Lm=0 


,a£ 


where, fi = y/( 1 — C 2 ); 


a = -C + tft 

* = VT~i); 
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Q‘2 + iP2 

and a m 
S m + iR m 


v Cm(a-a)+An^ 

„ho («- a ) 2+W m 

a + iuj m ] 

Dm 

a ^ + 2£a m + 1 


m 


Case (iii) /(t) = e at [E cos (fit) + F sin(Qt)] (iLo ^m^ m ) 

This type of non-homogeneous terms occur in equations (17) and (18) and it is 

similar to equation (19). Solution of equation (A.l) may be written as 


x{t) 

where, X\ (t) 
and x 2 (t) 

where, ft 

a 

and i 
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xi (t) + x 2 {t) 

e -Ct | (p 3 cos (/3t) + Qz sin(/3i)] / /5 
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for 2 < m < N. 
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for 0 < m < N. 
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Appendix B 


Laplace Transform of Equation 
(5.41) and Its Inverse 


The Laplace transform of Xq given in equation (5.41) may be written as 

A’oM = (B.l) 

where, N(s) = N 0 i(s)Yi(s)] iV ox (s) = (a x + 2( 1 s)(h n s 2 + h 12 s + h 13 ) ; 
h\\ = 1; h \ 2 = 4^(1 + o;)/E 2 ; h \ 3 — 2 a 2 (l + a) / R 2 
Yi(s) = 7 2 /(s(s + y) 2 ) for the rounded step; 7 = 7 S /T 0 ; 

D(s) = s 4 + 93s 3 + ^s 2 + <?is + ^oj 

93 = 2[(1 + a 2 + R 2 )(( 1 + C2) — 2a(Ci — C2)]/ R 2 ', 

0.2 — [1 ~ 2a;(a x — a 2 ) + a 2 + R 2 + I6C1C2)]/ R 2 ] 

Qi — 8[Ci a 2 + C2O1]/ R 2 ] Qo — 4aia2/ R 2 ■ 

The factorization of denominator D(s ) in equation (B.l) may be expressed as 
D{s ) = ipi(s)i> 2 (s)\ 

where, ipi{s) = [(s + a x ) 2 + $?]; lM s ) = K 5 + a 2 ) 2 + /?f]- 

Let us assume that A x = -cq + i/3 x ; A 2 = ~a 2 + i$ 2 \ are the complex roots of 
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the polynomial D(s) and Aj f- A 2 then the inverse Lapla.ee transform of Ad(.v) may be 


written as 


X 0 (T) = Ado + Adi (T) + X m (T) + X m {T) 
where, Ado = h n a } /qo\ 


A’oi (T) = exp(— 7 T)(A 0 i + Bq{T)\ 

B m = <3oio(-7); CoioM = -‘iow 51 ’ 

An = Qou(~ 7); Qou (•“>') “ Qoio(' S ) 

here, (') denotes differentiation with respect, to .s, 

Ad 2 (T) = exp(-Q !l T)(^02Cos(Ar) + 


B02 + i A2 = Qo 2 (Ai); 

let A 0 2o(Ai,A 2 ,T) = Ad 2 (T) then 


A 03 (T) = Ad 20 (A 2 ,A 1 ,T). 


The other parameters used in equations (5.39)-(5.41) an; 
A 14 — 1; h 15 = 4(i(l — a)/R 2 ] h i6 = 2a,(l - «)//?“; 
h 2 i = (1 - a)//? 2 ; h 22 = 4C ; 2 /^ 2 ; /i 23 = 2a 2 //? a ; 

A 2 4 = (1 + oi)/R 2 ] h 25 = 4Ci/i? 2 ; /i 2 c = 2a, /f? 2 . 
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